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SYNOPSIS 

STUDY OP PHOTOVOLTAIC DEVICES BASED OF 
METAL-SMIC OFDUC TOR STRUCTURES 

By 

BASABI BHAUMIK 

Department of Electrical Engineering 
Indian Institute of Technology, Kanpur, India 

August 1973 

Recently there has "been an intense interest in 
utilising photovoltaic conversion of solar energy for 
terrestrial applications. This has led to the exploration 
of novel structures which would provide an alternative to 
the conventional P-F junction solar cells. Our interest 
in this thesis has Been in studying some of these novel 
structures which are Based on metal-semiconductor junctions 
and act as Schottky Barriers or analogously related devices. 
Particular attention has Been paid to the understanding of 
the following aspects of these novel structures. 

1. The characteristics of the Schottky Barrier solar 
cells (SBSG) are rather well known [11. In the published 
literature, there are also some reports on the effect of 
temperature on the characteristics of these cells. Eor 
example, it has Been reported that the short circuit 
current Ig^ of silicon SBSO increases with increasing 
operating temperature. Following the concepts prevalent 
at that time (1975), the increase in I^q has Been attributed 
to the availability of extra photons, Belonging to the long- 
wavelength range of the solar spectrum, due to the lowering 
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of the forbidden energy gap of silicon with increasing 
temperature. Our study [2 ], presented in Chapter 1 of the 
thesis, however, shows that the contribution to the increase 
of Ig^ due to the mechanism mentioned above is negligible 
and cannot explain the experimental results. Instead, it 
has been found that good agreement with experimental results 
is obtained when the increase in with increasing 
temperature is calculated by taking the change of the 
absorption coefficient of the semiconductor with temperature 
at all wavelengths. This conclusion is not limited only to 
SBSCs but would be valid for P-N functions or other type of 
cells made on indirect bandgap semiconductors like silicon. 

2. The ma^'or disadvantage of SBSCs is due to their low 
open circuit voltage. There are two ways in which this can 
be obviated. One of these is by using Sehottky (MIS) cells 
[IJand the other is by using Sehottky (oppositely doped 
semiconducting interface) cells [3]. In the case of former 
cells, the performance critically depends on the thickness 
of the interfacial insulating layer. When this thickness 
is large, both the short circuit current and the open circuit 
voltage decrease leading to a decrease of the conversion 
efficiency. The mechanisms responsible for this behaviours 
have been found to be the variation of the ’illuminated 


barrier height’ and that of the density-of-surface states 
with the variation of the thickness of the interfacial layer. 
The details of this analysis are given in the second chapter. 



The other type of cells where an oppositely doped semiconducing 
interface layer is sandwitched between the metal and the bulk 
semiconductor have so far not received the attention they 
deserve. Keeping in view that the increase of effective 
barrier height, which is obtained due to the presence of a 
thin oppositely doped semiconducting interface, would be of 
interest not only from the point of view of solar cells but 
of many other devices where Schottky barriers are employed, 
we have developed a rather general method of analysis including 
the effect of mobile carriers for these devices [3 ]. This is 
given in the third chapter. 

3. All the devices mentioned above require a veiy thin 
film of metal to form the Schottky barrier because these films 
have to be transparent to light. This places stringent demand 
on the choice of the metal because it is required to simul- 
taneously meet the optical, electrical and mechanical criteria. 
A way out would be to use a lateral (gridded) structure where 
the metal film is not thin and continuous, but thick and very 
finely gridded. Since the carriers T ire photo generated in the 
area where there is no metal and arc laterally cbllected by 
the thick film of metal forming the Schottky barrier, we call 
this a lateral SBSG. Its analysis entails the solution of a 
two -dimensional continuity equation. A further complication 
arises because all the boundary conditions are not known. A 
method of analysis which assumes the functional form of the 
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boundary condition has been developed [4]and is given 
in the fourth chapter. It allows the determination of 
the boundary condition by proper matching and leads to the 
solution of the continuity equation. A brief summary of the 
properties of different types of cells considered in this 
thesis is given in the concluding fifth chapter. 
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CHAPTER 1 


SCHOTTKY BARRIER SOLAR CELLS 


Introduction 

The photovoltaic effect in semiconductor P-N junctions 
is well known and has been extensively utilized for power 
generation using solar energy in space. Lately there have 
been attempts to use this phenomenon for terrestrial generation 
of electrical power. The main impediment towards this, which 
again is well known, is the cost. To overcome this impediment 
two major approaches have been taken; first approach is towards 
cutting down the cost of P-H junction photovoltaic cells and 
the other approach is to search for novel structures which may 
prove viable. The work in this thesis is concerned with the 
latter approach I 1-5 J. 

The structures that have been investigated are based 
on the photovoltaic properties of Schottky barrier (metal- 
semiconductor) diodes and related devices like metal-(thin) 
insulator-semiconductor, metal-(thin)P-N, metal-(thin)N-P, 
metal-(thin) insulator- ( thin )P-N and metal-(thin) insulator- 
(thin)N-P diodes and are given in Ghapters 1, 2 and 5. A 
common disadvantage of all these devices is that they require 
a very thin metallic film which is almost transparent to light. 
This thin metal film has different lattice constant and 
different temperature coefficient of expansion from the 
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substrate and has to undergo temperature cycling during 
actual operation. Thus, there are chances that the metallic 
film -would peel off, affecting the durability of these types 
of cells. This can be obviated by evaporating a very closely 
spaced metalli c grid structure which can be thick, rather 
than evaporating the conventional continuous metallic layer 
which has to be thin [3,6,71. In the case of gridded 
structure, the carriers are photogenerated in the exposed 
region of the semiconductor and are laterally collected by 
the metallic strip forming a contact based on any of the 
structures mentioned above. Such a cell has been called a 
lateral cell and, as is obvious, it entails a two-dimensional 
analysis [33. This is the subject matter of Chapter 4. 1,J e 
start the considerations of this chapter with a study of 
Schottky barrier solar cell. 

1 . 1 Schottky Barrier Solar Cell (SBSC) 

A good review of Schottky barrier solar cell (S3SC) 
is available in the book by Hovel [8l . Here we mention tnose 
aspects which are either required later or are necessary for 
the sake of completeness. A particular aspect needed later 
concerns the temperature variation of the characteristics of 
solar cells. Though the development given in this chapter 
is in the context of SBSC, the arguments are general enough 
to hold in the case of other photovoltaic structures also. 
The consideration of dependence of performance on temperature 
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has been necessitated by the fact that in actual operation 
the operating temperature of a solar cell keeps on changing. 

A perusal of the literature showed that there were some 
ambiguities regarding the physical mechanisms responsible, 
for the dependence of the : characteri sties of SBSO on 
temperature. A clarification of these II ] is given in 
Section 1.3, and is preceeded by a brief historical comment 
and the derivation of' the current voltage characteristics 
of - SBSC . 

The photovoltaic properties of those metal semi- 
conductor contacts which form a Schottky barrier have been 
known for a fairly long time [ 9 ) . At an earlier stage of 
development these were mostly used as photodiodes, fhe 
basic operative mechanism of this system as photovoltaic 
device for energy conversion has been discussed by fehovec 
[10 ] in 1943. However, barring these sporadic efforts, till 
recently no systematic attempt had been made to correlate the 
output characteristics of these devices with the process of 
fabrication. The recent resurgence of interest in the 
photovoltaic conversion of solar energy and ensuing search 
for photovoltaic structures which would provide alternatives 
to conventional P-F junction solar cells have, however, changed 
this picture. An important step in this direction was the 
theoretical analysis of Pulfrey and McQuat till which showed 
that under ideal conditions the efficiency of a Schottky 
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barrier solar cell' (SBSG) would be comparable to that of a 
P-N junction solar .cell-. Added to this is the advantage 
that the fabrication of SBSC would be simpler and processing 
could be done at lov;er temperature. In actual practice, 
however, the conventional SBSC has not lived upto the promise 
due to its low open circuit voltage and the ensuing lower 
efficiency. The structure that seem to be viable are based 
on the analogous metal-(thin) insulator-semiconductor and 
are called tunnelling MIS (or MOS) diodes. These have higher 
open circuit voltage and efficiency as compared to SBSC. We 
would take up their study in the next chapter after studying 
the behaviour of SBSC, which is simpler, in this chapter. 

A convenient, but round-about, way of coming to the 
operative mechanisms of 33SC is through those of a P-N 
junction solar cell which are more widely known. The energy 
band diagram of a P-N junction with light incident on it from 
the left is shown in Pig. 1.1(a). Here, the photons which 
have energy larger than the bandgap of the semiconductor create 
electron-hole pairs. The electric field in the region of the 
metallurgical junction acts to separate the pairs which are 
generated in the depletion region, and provide drift field for 
the minority carriers which manage to diffuse to the edges of 
the depletion regions on both P and N sides. Finally the 
separated carriers are collected by the two ohmic contacts, 
one on the back side (N side) being continuous and the other 
on the front (P side) being gridded to let the light in. In 
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summary, the operation of a photovoltaic cell requires the 
photo -generation of electron-hole pairs, their separation 
by an electric field and the collection of these separated 
carriers by the metallic contacts so that they can flow, 
through the external circuit. In addition one of the 
metallic contacts should be in position to let the light in, 
i.e. either it has to be very thin or of gridded structure. 

For a moment let us now imagin^that the P-EF junction 
is sectioned at plane AA* and light is incident from the 
left of plane AA'. When one considers the possibility of 
using this structure for photovoltaic conversion one finds 
that it meets the necessary conditions of generation of 
electron-hole pairs and their separation by an electric 
field, but does not meet the requirement of proper collection 
of the photogenerated carriers. Of course, there is the 
question that in the absence of P semiconductor how has the 
band bending to the right of AA’ been obtained. Both 

this and the requirement of collection can be simultaneously 
met if one puts metallic contact at AA’ so that it forms a 
Schottky barrier as shown in Pig. 1.1(b). It is needless to 
say that this metallic contact should be almost transparent 
to the incident radiation,. Phis is the basic principle of 
operation of a Schottky barrier solar cell (SBSC) based on 
metal-N type structure. If we repeat. our argument by looking 
to the left of AA' and by considering light incident from 
right, a metal-P type S3SG can be observed. 
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The considerations given above clearly show the analogy 
between SBSG and P-N junction solar cells. Thus they allow us 
to carry through, but with caution, many of the concepts 
developed for the latter in the case of SBSG also. 

1 .2 I-V Characteristics of SBSG 

Here there are two cases of interest, (i) The diode 
is forward biased by an applied voltage Y and an I-Y relation 
ensues . There is no illumination hence we would call this as 
the characteristics of the 'dark' diode, (ii) The diode gets 
forward biased at a voltage V by the connection of a load 
in the presence of illumination. Part of the I-Y characteristic 
is in the fourth quadrant and would be denoted as the characteris- 
tics of the 1 illuminated' diode. We take up the case of dark 

f 

diode first. 

1 .2.1 Park current of Schottky barrier 

Current voltage characteristics of dark Schottky 
barriers have been well studied in the published literature 
I 123. Most of these studies are limited to small area diodes 
so that the fluctuations of any parameter on the diode do not 
complicate the properties. By necessity, solar cells have 
to be of comparatively larger area and due to this some new 
considerations come in. These considerations have been 
briefly pointed ou t after first giving a method to determine 
the mechanism of current transport in diodes which have 
homogeneous properties. 




It is by now well established that depending upon the 
barrier height, doping, temperature and applied bias the 
mechanism of current transport through a Schottky barrier 
diode can either be thermionic emission (IE), thermionic 
field emission (TEE) or field emission (EE). It would be 
of interest to find a method so that by making measurements 
on a finished Schottky barrier one could determine as to 
which one of the three mechanisms mentioned above is operativ 
in the diode. Such a method has been suggested by Saxena C 131* 
The method consists of experimentally plotting a graph of 
nT versus I where n is the ideality factor of the diode, A 
typical plot of this is shown in Eig. 1.2. The graph I 
corresponds to the case of an ideal thermionic emission 
with n=1 , whereas graph II also corresponds to thermionic 
emission but with n > 1 . There is one more case shown by 
graph III where thermionic emission is the mechanism, but this 
graph does not pass through origin, rather id? is parallel to 
n=1 line. This is the Case of 1 T f anamoly and has been 
discussed by several workers [14-9* 14b]. 

In the case that the thermionic field emission is 
the mechanism of current transport through the diode under 
consideration one obtains graph IV as the plot, whereas in 
the case of the field emission graph V results. Thus, one 
sees that the plot of nf versus T provides a simple method 
for knowing the mechanism of current transport. However, the 
perusal of recent literature shows that there have been some 
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Fig. 1.2 Plot of V versus (KT/q) to determine the mechanism of current transport 
in Schottky Barriers (See Text) 
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ambiguities in this regard. A case in point is the experi- 
mental results of Vernon and Anderson [ 15 J. They have fabri- 
cated large area Schottky barrier solar cells and in order to 
determine the mechanism of current transport in their samples, 
they have plotted nT versus 1 from the-I-V characteristics of 
SBSO. They have found this plot to be independent of tempera- 
ture and have concluded therefrom that field emission is the 
dominant mechanism of current transport. This conclusion, 
however, seems to be suspect because neither is the semicon- 
ductor degenerately doped nor is the temperature low enough 
for field emission to dominate in current transport. It is 
true that if one only considers the conventional mechanisms 
(thermionic emission, thermionic field emission, and field 
emission) of current transport in a Schottky barrier diode, 
one is left with no choice but to conclude that the mechanism 
in this case is field emission. It should be noted, however, 
that the SBSO’s are large area devices with no guard rings. 

In a large area device it is quite possible that the barrier 
height is not uniform over the whole sample. Johnson etal.f 16?, 
in a different context, have investigated the effect of the 
variation of barrier height on the I-V characteristic of a 

metal-semiconductor junction. In some cases they have observed 

tT\V «.?>$<£ djL 

that the plot of the product of the temperature and the slope 
of the I-V characteristic versus temperature is independent 

A 

of temperature and have ascribed this to the variation of 
barrier height. Although the argument is not conclusive, it 
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may be that this rather than field emission is the controlling 
mechanism of current transport in large area SBSO’s. 

Thus it is evident that in the case of large area 
Schottky barrier diodes great care must be exercised in using 
the method of Saxena to determine the mechanism of current 
transport. This has been followed through by another worker 
of our group [ 171. In this thesis, however, unless otherwise 
specified, we would take the mechanism of current transport 
to be thermionic emission and use the equation: 

I = I [ exp(qV/nkT) -1] (1.1) 

where I = aJ , a being the area of the diode and 
o o 

** 2 

J 0 = A t [exp(-q0 B /kT)] (1.1a) 

for the dark current of SBSO. Here A** is the Richardson 
constant, 0 B is the barrier height, T is temperature, q is 
the charge of electron and k is the Boltzmann constant. 

The value of n would be taken to be unity unless otherwise 
specified. 

1.2.2 I-V characteristics of S3SC 

The typical nature of current voltage characteristic 
in the presence of illumination is shown in Big. 1.1(c). At 
a given load R^, a voltage Y is developed across the cell. 

For the present let us assume that the effect of the series 
and shunt resistances and the recombination of photo -generated 
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carriers in the depletion region and at the interface are 
negligible. Under these conditions, the following analysis 
holds. At the voltage 7, let the current be I and the 
current density be’ J for a given area a. Let us consider a 
metal-1' semiconductor Schottky barrier cell. At any cross- 
section of the cell, we have 

J = Jp(x) + J Q (x) (1.2) 

where Jp(x) is due to the flow of majority carriers and 
under the assumptions given above, can be obtained from 
eqn.(l.l). J (x) due the flow of minority carriers 
and is determined below. In general, it is given as: 

J n (x) = J qD (x) + J nQ ( x ) + J R1 ( X ) + j r 2 (x) (1.2a) 

where J n p( x ) is the current due to the photo-generated 
carriers in the depletion region, J^q(x) is due to diffusion 
of photo-generated carriers in the quasi -neutral region just 
beyond the depletion region, Jj^( x ) an< ^ ^ are ^ ue 
recombination of carriers in the depletion region and at the 
interface, respectively. As mentioned earlier, here we would 
neglect the contribution of J R ^(x) and J-^* '■® ae ex P res8 i° nB 

for J n q(x) and J^p(x)are developed below. Under ©commonly 
employed assumptions [ IS ] J n ^(x) can be obtained by using 
the continuity equation for the excess minority carriers 
Up which is given as (1 ] 
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d n 

® ^ + Q(X) {1-R( \1)a(X )exp { -a(X)(x+b))= 0 (1.3) 

n dx 2 T n 

Here b is the width of the depletion region (Pig. 1.3a) and 


other symbols are explained in th~ 


First 


of all, let us consider eqn.(1.3) when the cell is short 


circuited. 


The boundary conditions are 

at x = b; Up = N 
dn 

at x - W; D n d^i x=W = " s °p< W) 


(1.4) 

(1.5) 


Here s is the surface recombination velocity at x=W and W 
is the total width of the semiconductor. i. comment on the 
choice of the boundary conditions at x = b is given later. 
For the present we note that the solution of eqn.(1.3) with 
constraints of eqns. (1.4) and (1.5) is given as 

n = c 1 exp(x/L n )+ c 2 exp(-x/L n )+ 9 (1.6) 


- a 2 ) 


where 


c^ = X"C t (a D n ~s)exp{ -(^- + W) } + ( s + ^p) exp(-ab) 

n n 

exp( ^ 2 - ) 3 / 2A - exp(-( +ab> )} + N exp(- ) 


{ 1- t (s + yS.)/2A] exp(^f^) ) 


(1.6a) 


: 2 = ~X ^ ( a D n - s)exp(- ~ ) + (s+ ^S-) ex P(- a "b+ /2A 


+ N(s+ T S)exp(W/L )/2A 


(1.6b) 


where 

X = Q(1-R) a / - O. 2 ) 


(1.6c) 



LU3\ DM 



Fig.1.3(a) Energy band diagram of a P- semiconductor diode showing the quasi fermi levels, 
(b) Plot of E Fn (b) and N versus J n (or N a ). (Seetext). 
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D 

A = s sinh [ (W-b)/L ] + cosh. [ (W-b)/L ] (1.6 ct) 

-u n n 


Prom eqn. (1.6) one obtains the short circuit current IgQ^( x) 
as: 


I SGQ 


(X) 


dn 

qaL 

H n dx 


x=»b 


X> 

aL exp( -a b) (a L -s) y - aL 

= qaQ(l-R) 3 - [ + A M } 

(1- a L n ) A exp (a W) exp (ab) 

(1.7) 

2 

where L = D t , a is the area of the cell and 
n n n 


+0 


B = s coshf (W-b)/I n 1 + (D/V sinh[ ^ ¥ "' b ) /L n 3 ^* 7a ) 


G = 


qaD n N 


[1 


D 

_£L( 

A v 3D 


S In/ W_;bs , 

+ y— ) exp( — T - ) ] 


(1.7b) 


n n n n 

Note that inadvertently a printing error in eqn. (1.7b) has 
appeared in Ref. til . The corrected equation is given here. 


Before proceeding further we \vould like to comment on 
the component 0 of short circuit current. It is customary in 
the case of short circuited cell to take N = 0 in the boundary 
condition at x=b (eqn. 1.4), in which case 0 becomes zero. The 
fact that N is non-zero even in the case of a short circuited 
cell can be seen as follows [19-21] . Let us consider the 
energy band diagram of a metal P semiconductor Schottky barrier 
solar cell shown in Pig. 1.3(a). The considerations for the 
metal N-semiconductor SBSG are very similar. At the metal- 
semiconductor interface the minority carrier _imref IL 
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is shown to merge with the metal Fermi level Ep^, indicating 
an infinite recombination velocity. Poliowing an analysis 
similar to that of card 1 1 9 ^ , it can be shown that the 
minority carrier imref Ep^ does not remain constant in the 
depletion region as is commonly assumed. Let us take the 
reference of energy at the metal Fermi level and denote the 
value of E^ at the edge of the depletion region at x = b by 
E Pn(b)‘ minori ' t y carrier current density in the absence 

of recombination is given by 


d E. 


J = q n V 
n H n 


Fn 


dx 


(1.8a) 


where 

E r -E 

n = exp(- -T_ ™) 
0 kT 


The symbols are defined in Appendix A. 


( 1 . Sb ) 


Using eqn. (1.8), rearranging, and integrating for 
0 < x < b one gets 


b J Fn v u 

/ exp(E r /kT) dx = j exp(E- /kT) dE (1.9a) 

o N 0 ^n 0 E-. (o) En Pn 

-Til 


Evaluating this, and making some simplifying assumptions, 
one obtains : 


expO 


K _ 0 >) - 


Fn'' 1 Fn V \ 4 T ... Jh 


kT 


)-1 = J n #(1/ n i u n ) (ire g /2q kT)^ exp[ 

(1.9b) 


-E^(o) 


x^n 


TF 


•3 


In metal P semiconductor Schottky barrier diodes, Ep n (o) = 0 


Using this, one obtains the equivalent expression of Card t19l 
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In the case of metal F-type contacts, eqn. (1.9b) still holds 
with and y n replaced by -B ^ , F^ and y p respec- 

tively . 

A plot of Ep n (b) as a function of photo-generated 
current density J n using eqn. (1.9b) is given in Fig. 1.3. 

Also shown on this graph is the value of II given as 

F = Fq exp t-(E Q -E ;pn (b) - qV p )/kT] (1.10a) 

= a po exp [ E Fn (b)/kT 1 (1.10b) 

where qV p is the difference between Eermi level and the 
valence band. The existence of non-zero value of F in the 
case of short circuited photodiode was first noted by 
Berry [20] in the case of P-F junctions, and the numerical 
value of F was deduced from the short-circuited photo- 
capacitance of commercial solar cells. The value obtained 

1 4 -3 

in this case was F = 10 cm . Moore [21] suggested that 

the analytical model of Berry based on depletion region 

photo-capacitance was not correct, rather a diffusion model of 

photocapacitance should be used. This correction gave a 
11 -3 

value of F - 10 cm . The point to note is that if one 

11 -3 

evaluates the current component 0 using N - 10 cm then 
the contribution of the term is negligible compared to the 
other terms. Hence the presence of F in the case of 
Schottky barrier solar cells is only of academic interest. 
However, this term becomes of considerable significance in the 
case of MIS solar cells as would be shown in the next chapter. 



IS 


There we would use an expression for the current 
density which is equal to (-C/a) where 0 is given "by 

eqn.(1.7b). If we evaluate the current C by assuming that 
¥ >> I' n such that the recombination velocity s at the 
contact x = W does not come into picture, one obtains a 
modified form of eqn.(1.7b) as 


q a P 

0 - sy 2 ' < s - n P o ) 


Use of eqn.( 1.10b) in this expression gives an expression 
for > which has been defined above, as: 


q -D 


n 


J diff = " G/a = u po 1 r 


{ exp(- 


W*>- 


-) -1 } 


(1 ,11a) 


n 


In the case of a metal ft semiconductor Schottky barrier 
diode eqn.(l.lla) gets modified as: 

P E (b) 

J diff = ~°/ a = 4 T 2 { ex P < &F— > - 1 } (I.Hb) 


L p ‘no 


kT 


where E (b) < 0, if metal Fermi level is taken to be the 

j? p 

reference energy. As already mentioned, these equations 
would be used in the next chapter while analysing Schottky 
(MIS) solar cells. 


The current generated in the space charge region 
(0 < x < b) when recombination is neglected is given by 


I 


SOP 


( X ) 


q a { 1-R( X ) } 


b 

J Q( X ) a ( X ) exp{- a( X )x } dx 

(1.12) 


The total current Ig^(X) for a monochromatic light of 
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wavelength A is given bytthe sum of eqns. (1.7) and (1.12) as 

I so ( A ) = I 30Q ( A ) + I gCD ( A ) = a J L ( A ) + 0 (1.13) 

This current has "been evaluated for a metal I -silicon 
Schottky barrier solar cell and the values are given when 
we calculate the efficiency of the cell. 


Having obtained the short circuit current let us 
obtain the open circuit voltage. This can be obtained by 
equating the sum of the majority and minority carrier currents 
to zero and solving for the voltage. The majority carrier 
component can be obtained either by using the thermionic 
emission model or by diffusion model of current transport. 

If thermionic emission model is used then the use of eqn.(l.l) 
and (1.3) gives the total current as 


I = I Q [ exp(qV/kT) -1 } - I gc 


(1.14) 


3y equating 1=0 and solving for V one gets the well known 


expression for Y^ as 


n kl 


V 


00 


In [ 1+ -JO 3 
*0 


(1*15) 


An expression for Y^ using the diffusion model of 
the majority carrier transport has been obtained by Dubey 
and Paranjape [ 22]. They have used a method of boundary 
matching and obtained the expression of Y^ in terms of a 
recombination velocity at the metal semiconductor interface. 
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However, for the parameters that are com.ionly useful in 
solar cells, it is more likely that thermionic emission 
model rather than diffusion model would hold true, hence we 
would use eqn, (1.15) for V QC . 


To complete tne analysis, we want to obtain an 
expression for the efficiency of the cell. This can he done 
in the conventional manner using eqn. (1.14). The power output 
is given by P = IV and the voltage V , a t which the maximum 
power is available, is obtained by solving for 3 P/37 = 0. 

This leads to a transcendental equation 
nV qV I 

(1 + -#) expt-gSE) _ (1 + Jj£) 0.16) 

where Iq is the reverse saturation current of the diode and 
IgQ is obtained by integrating eqn. (1.15) over the range of 
wavelength of interest. 


Using the value of V obtained from eqn. (1.16) in 

eqn. (1.14) the value of the current I at the maximum 

max 

power point can be obtained. The expression for the effi- 
ciency ( r) ) is now obtained as 


IV IV 

„ max max r-n -a \ SC OC 

n = = (p.p.)— p— 


(1.17) 


' in 


m 


where is the intensity per unit area of the incident 
light and P.P. is the fill factor, which is given by the 
ratio 

3r.r. = \ax (1.18) 
i sc oc 
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It is important to point out here that the theoretically 
calculated efficiency of the same system using eqn.(1.17) can 
vary "between wide limits depending upon the assumptions made and 
the values of parameter assumed. For example, for cells made 
on IT-type silicon the values of efficiency ranging from 
approximately 20 $ til] (obtained for the equilibrium barrier 
height 0-g of the order of the energy gap Eg and assuming 
negligible recombination and unit quantum efficiency) to 
5$ {23 J (in the case of gold silicon contact taking incomplete 
absorption and some recombination into account) has been 
reported. Thus, caution must be exercised in completiy 
specifying the conditions under which a number for efficiency 
has been obtained from the theoretical calculations. It is 
also worthwhile to note here that during the earlier stages 
of the development of SBSC’s, one very intriguing feature 
had appeared. It had been found that the experimentally 
obtained efficiencies where higher than theoretical efficiencies 
calculated taking realistic assumptions into account. It was 
soon realized that this situation had arisen because the 
experimental cells were not actually Schottky barrier cells 
for which the theoretical model was built. The presence of 
the thin oxide layer formed during processing, between the 
metal and silicon^ had given rise to tunneling MIS solar cells 
which requires different modelling than S3SC’s. A calculation 
of the efficiency for a SBSO made on P-type silicon and 
operating at 300°K has been done as follows: 



22 


(i) Ig q( X) has been evaluated from eqn. (1.1 3) using 

2 

the following values of parameters, D n = 10 cm / sec., 


L = TOO s = 10 cm/ sec* f W = 200 ynij 


1 cm, b = 1 pm. 


Q( X ) corresponding to AMO sunlight [ 24 ] and the values of 
a( X) from Dash and Newman [25] have been used. The total 
short circuit current has been obtained by integrating Ig£,(x) 

from x = 0.4/4111 to x = 1.1 /4m. This calculation gives 

2 2 
IgQ = 25 mA/cm f 0% R = 0.1 and 16.5 mA/cm for R = 0.4f 


(ii) Before calculating Y nn from eqn. (1.15), we want to - 

(JO 

note some of its features. If (Ig^/Ip) » 1, and Iq = aJ^, 
where J Q is given by eqn. (1.1a), 


Y oc = n0 


nkT r 
■B + — ln [ 


■SC 


aA** T‘ 


] 


n 


2 2 ** 

Taking Ig 0 = 25 mA/cm , a=1 cm , A 
= 1 , and T = 300°K, one gets: 


(1.19) 


32 Amp/cm^/°K^, 


V 00 = " °*5 volts 


If one takes 0^ = Eg = 1.12 eV, then Y^ = 0.62 7 and if 
0g = O.S eY then = 0.3 "V. Thus it is apparant that 
for realistic values of barrier heights the open circuit 
voltage of SBSC would be very small as compared to that of 
P-N junction solar cells. This is a serious limitation of 
SBSC and has limited its usefulness. 

(iii) The calculated values of the efficiency of SBSC would 
depend upon the choice of 0g. If we take 0^ = we get 
n = 6. 5% for R = 0^0 which is an upper limit on the 
efficiency in this case. 
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1 . 3 Variation of Parameters of S3SG with Temperature 

The main reason for the interest in the study of 
temperature dependence of performance of a solar cell is 
due to the fact that in actual operation it has to undergo 
temperature cycling. Naturally this problem had attracted 
attention at an early stage of the development of solar cells 
[26]. Recently, Vernen and Anderson [ 15] have experimentally 
measured the temperature dependence of the short circuit 
current, open circuit voltage and some other parameters of 
silicon SBSO' s. Let us consider first the variation of the 
short circuit current IgQ* It has been found that it increases 
with increasing temperature [15]. The explanation given for 
this increase is that as temperature increases the bandgap 
Eg decreases, thus making the absorption of some more photons 
of longer wavelengths, which were unutilized earlier, possible. 
Particularly for their data, Vernan and Anderson have suggested 
that "the increase in Igg can be accounted for by 100°0 temper- 
ature increase producing a 0.03 eV decreased energy gap and 
0.03 um increased useful wavelength range". That this is not 
the case can be physically seen from the following arguments. 


T 1 

Let be the wavelength of the photons corresponding 

T 

to Eg at temperature T^ and ^ be the wavelength of the 
photons corresponding to Eg at temperature Tg- 

T 1 

Eor silocon: =1.1 Pi at 300°K 

*2 = 1.13 hm at 400°K. 


and 
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Here use has been made of the expression [ 27 ] 

2 

a T 

E p (T) = E p (0) — (1.20) 

G G (T + 8 ) 

where for silicon a= 7.02 x 10 ^(eY/°K); $= 110S°K 

and Eg(0) = 1.16 eV. This expression for Eg(T) is valid for 0° 

to 900°K. 


The value of the absorption coefficient at A- 1.1 ym 
for silicon is around 10 cm [25], which means that the 
photons at these wavelengths would be absorbed at a depth 
of 1000y m. For good silicon crystal the diffusion length is 
of the order of 100 U m. Hence it is evident that photo- 
Carriers generated at these wavelengths can not contribute to 
the short circa t current and thus do not account for the 
experimentally obtained value of dlgg/dT. Instead, in the 
present work, it has been found til that for the range of 
wavelengths ( X) of solar radiation useful for photogeneration 
in silicon, the increase in the absorption coefficient a with 
increase in temperature is the dominant mechanism which gives 
rise to the increase of Igg with temperature. To show this 
the expression for Igg( X ) as the sum of eqns. (1.7) and (1.11) 
is used. At each temperature, I g( X) is integrated over the 


whole useful wavelength range. By repeating the process at 
different temperatures, the temperature dependence of Igg is 
obtained. The results of the computations are shown in Fig. 1.4. 



Fig. 1.4 Xgg v«rsu« T,wifeh L n as pcra»«t#r 
Xnast shows scharaatio of 8B8C. 
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The computations for Ig^( A ) have been done for an 
2 

area a of 1 cm and between the temperatures 300° and 425 °K. 

_ 2—14 

D„ and S have been taken r23i to be 10 cm sec - and 10 cm 
n n 1 1 

sec - ' 1 , respectively. L n is treated as a parameter ranging from 
10 to 100 n m. The width W = 200 pm and the width of the 
depletion region b = 1 V m. The temperature variation of E^ 
for silicon is already given by eqn.(1.20). Q( X ) corresponding 
to AMO sunlight has been used, and it has been assumed that the 
transmittance of the thin metal film is unity over the whole 
range and there is 40$ reflection at the metal semiconductor 
interface. The value of the absorption coefficient a (X ) at 
different temperatures has been evaluated by using the equ ition 


t-29, 30 ] 


a = 


(hv - E n + K 8. ) 2 

E [ £ 

i expCe^T) - 1 


(hv - E - E e.) 2 

+ — ] 0 . 

1 - exp(- 6^/T) 


(1.2l) 


Here hv is the energy of photons, K0 ^ is the energy of 

phonons, i is the index for different types of phonons 

participating in the process, andG-^is a constant. We have 

3 -1 -2 

taken bh = 4.5 x 10 cm eV over the whole range 
wavelengths and temperatures. This value of 0^ has been 
obtained by comparing eqn.(1.21) with the data of Dash and 
Newmann [25] and by considering the proper threshold energies 
for different type of phonons as given by He lean [311. The 
value of in the range 1.175-1.21 eV has been taken to be 
212°K corresponding to transverse accoustical phonons, and 
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in the range 1.21 - 3.0 eV to be 670°K corresponding to 
transverse optical phonons. Thus the absorption coefficient a 
has been computed in the range 300-425° I, "’hich are also the 
respective ranges in the calculation of 


The calculated plot of I,,- versus temperature is 
shown in Pig. 1.4. The following features should be noted. 

(a) The dashed curves in Pig* 1.4 do not show any variation 
with temperature. They have been calculated by changing the 
integration limit to take care of the lowering of E^ with 
increase in temperature, but the value of a ( X ) is kept 
fix^d to the value at 300°iC. Thus, it is not true that 
"the increase in can be accounted for by 100°0 temperature 
increase producing 0.03 eV decreased energy gap and 0.03 y m 
increased useful wavelength range. This can be seen physically 
in the following way. The absorption coefficient in silicon 
for photons corresponding to E^ at 300°I is 'v 10 cm . T, ftien 
Eg decreases with increasing temperature, the number of excess 
photons absorbed is negligible due to tie extremely low value 


of a in this range. 

(b) The solid lines in Pig. 1.4 show variation with 

_2 

temperature and dI S( g/dT ranges from 4 to 7 y A cm /°K for 
values of l n from 10 to 100 u m. In these calculations 
temperature variation of a , as described earlier, has been 


incorporated along with the proper change in the limits of 
integration to take care of the variation of E^ with tempera- 
ture. The values of dl^/dl in Pig. 1.4 are comparabl e with 
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2 

the experimental value [15 J of 5 yA cm” /°X. An exact 
agreement could only be sought if the parameters of the 
experimental samples were known exactly. It is to be noted 
that the considerations given here are not limited only to 
S3SG. We consider this to be an important contribution of 
the present work. 


Having studied the temperature dependence of Ig^, it 
would also be of interest to know tne temperature dependence 
of the open circuit voltage which is given by eqn.( 1 . 19 ). 

We observe that the evaluation of this equation at different 
temperatures requires a knowledge of the variation of 0-g and 
n with temperature. A survey of the published literature has 
shown that (d0T-,/dT) and (dn/dT) are sensitive to the method 

JD 

of fabrication and this has given rise to several conflicting 
reports. This situation has been succinctly summarized 
by Hackam and Haroop [ 32]and we reproduce some of their 
comments in Table 1.1. A perusal of this shows that d0^/dT 
and dn/dT should be found for a particular device and using 
these the value of (dV^/dT) should be calculated* The same 
uncertainties that affect the calculation of dV Q( -,/dT also 
affect the uncertainties in (dn/dT), , hence in these cases 
we have not carried out the computations any further. 

To conclude, in this chapter, we have obtained the 


expressions for the current voltage characteristics of 
SBSCs’ and have studied some aspects of their temperature 
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Table 

1.1: Temperature dependence of 0 B and n in 

silicon. 

Metal-Si 

systems 

d0g/dT and dn/dT 

References 

Gr-nSi 

Same values of 0,> and n at 77°K, 196°A 

i 

and 289°X are obtained by — n— -V method 

Gowley etal, 
IEEE Trans. , 
ED-15,761,1968 

Titanium-Si 

(both P and 
N-type) 

The barrier heights were determined from 
measurements of the reverse current 
density and activation of reverse c 

current. d0-/dT varies from 0 to -IxlO - ' 5 
ev/°£ near 300°K 

Cowley, Solid- 
St. -Electron. , 
12, 403, 1970. 

Au,Ag and 

Gu on n-Si 

d0 B /dT = -(2.67+0.60) x 10"* 4 eV °K 
and n varies from 1.01 to 1.02. 

Arizumi etal., 
Jap-.J.Appl. 
Rhys, 8, 749, 197; 

Au-nSi 

Equality of temperature dependence of 

0 B and E & 

Crowell etal. ; 
App. Phys -Lett, 
4,91,1964 

Au-nSi 

d0g/dT 'v o from capacitively determined 

0^; d0g/dT n, large for 0^ determined 
by photoelectric technique 

Xahng, Solid 
s t. Electron. , 
6,281 ,1963 

NiGr-nSi 

n depends on forward bias voltage. At 
low voltage, n decreased with increasing 
temperature while at higher bias n was 
independent of temperature 

Yu and Snow, 
J.Appl.Phys. , : 
39,3008,1968 

Gr-nSi 

Both 0 B and n depend on bias voltage 

Levine, J.Appl 

Phys. ,42,3991 . 
1971 

Ni-Si 

Au-Si 

(i) d0g/dT >0; C-V method used. 
Temperature range used is 77°X- 

330Oi - a. 

(ii) d0 Q /dI o from active energy 

Jj k 

Padovani, J.Ap, : 
Phys., 33,391 ■ 

1967. 


plots. But this plot yields higher 
value of Richardson constant. When 
d0g/dT = -2 x10-4/o:{ is invoiced , 
accepted literature value of 
Richardson constant is obtained. 
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dependence. As has already been mentioned, the SBSGs’ 
have not lived upto their promise so far efficiency of 
energy conversion is concerned .This is mainly because of 
their low open circuit voltage. Tunneling MIS solar cells 
have shown better performance and are the subject matter of 
the next chapter. 



CHAPTER 2 


SCHOTTKY (MIS) SOME CELLS 


2. 1 Introduction 

It has "been found both experimentally [1-5Jand 
theoretically [6-15] that interposition of a thin insulating 
layer between the metal and the semiconductor of a Schottky 
barrier solar cell increases the open circuit voltage, thereby 
the efficiency. This increase, however, critically depends on 


the thickness of the insulating layer, <5 . A typical experi- 
mental curve, obtained by Lillington and Townsend [ 1 ] exhibiting 
the dependence of the open circuit voltage and short circuit 
current Ig^ on 5 , in the case of a Au-SiOg-nSi cell, is given 
in Pig. 2.1. The notable features of this plot are the following: 
(i) Igg remains constant with 6 upto a critical thickness 


SgQ and then decreases. 

(ii) Vqq increases with 6 , reaches a maximum and flattens 
out and then at a critical thickness 6 ^ starts 


decreasing. 


Although some heuiristic arguon nts have been provided 
in the published literature to explain the behaviour of Pig.2.1, 
to the best of our knowledge no analytical attempt has been 
made to quantitatively obtain the behaviour of Vqq f° r 
6 > 6 and 1^ for 6 >< %q* In this chapter we attempt 

to develop such an analytical theory. This has been possible 
because of the inclusion of (a) The variation of the density 
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of surface states with 6, while calculating the interface 
recombination current and the barrier height and (b) The 
shift, under illumination, of the minority carrier quasi- 
Fermi level from the metal Fermi level. 

In the published literature there are several [151 
theoretical analysis of MIS solar cells. Each of these has 
emphasized different aspects of the problem. In particular, 
we would like to mention here th e works of Landsb-erg and 
Klimpke [12], Buxo etal.[13], Card and Yang[11] and Card [14]. 
Landsberg and KLimpke have considered the effect of recombina- 
tion at the interface in detail and have investigated the effect 
of the density of surface states on the efficiency of the cell, 
but they have ignored the splitting of the minority carrier 
quasi-Permi level from the metal Fermi level. This latter 
effect has been taken into account by 3uxo etal.[13 J, but 
these authors have ignored any recombination at the interface 
or any other effect due to the density of surface states. Card 
and Yang [11] were perhaps the first to point out the role 
played by the splitting of the quasi Fermi level of the minority 
carriers from the metal Fermi level, but they have restricted 
themselves to a qualitative [11 ] or a semi-quantitative [14 ] 
analysis. The work of some other workers in the area has 
already been reviewed by Olsen in a recent paper [15] hence 
is not being repeated here. In the present work attention has 
been restricted mainly to the analysis of the variation of the 
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short circuit current and the open circuit voltage with the 
oxide thickness. To do this both the factors (a) and (b) 
mentioned above have been included. Though several workers 
have considered them separately, their joint effect, to the 
best of our knowledge has not been investigated earlier. It 
is also to be pointed out that the importance of variation 
of the density of surface states with the oxide thickness in 
explaining the behaviour of MIS cells has not been realized 
in any of the earlier published works. 


The expressions for Y^q and Ig~, with the inclusion 
of the factors (a) and (b) mentioned above have been develop- 
ed in Sec. 2. 3, where it has also been shown that with 
suitable approximation in the expression given here, the 
familiar expressions for given in the literature by 

other workers fall out. In order to obtain numerical values 
from these expressions one would need the value of the barrier 
height under illumination and the values of some other para- 
meters. The expressions for these have been developed in 
Sec. 2.2. The chapter has been concluded by developing an 
iterative scheme, which though complicated, allows the 
calculations of the illuminated current voltage characteris- 
tics of the Schottky (MIS) solar cell. 


2.2 Energy Band Diagra m 

Our main interest is in developing the expressions 
for the open circuit voltage and the short circuit current 
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of an MIS cell. Before doing it, it i s necessary to discuss 
the features of the energy band diagram of an MIS solar cell. 
These band diagrams (i) in thermal equilibrium, (ii) under 
illumination (such that a forward voltage develops when a 
load is connected) and (iii) under forward bias (but without 
illumination), are given in Fig. 2.2. The semiconductor is 
taken to be N-type but similar considerations would also hold 
for F-type semiconductor. In the following we consider the 
diagrams of Fig. 2. 2 one by one. 


2.2.1 Thermal Equilibrium 

In thermal equilibrium, the interest is in knowing 
the barrier height 0^. This has been calculated by Cowley and 
Sze [16] taking the occupation probability as a step function, 
i.e. a 0°K approximation; and by landsberg and Klimpke [12 ] 
by taking the occupation probability as Fermi Dirac distribu- 
tion function and is given below (see Appendix 3, equation 
( 3 . 16 )) 


0 B * f0M~X> 


M 


S 2 P 


2 e i 


q 2 6 

3 0 

e s’o 


q V 


1 +exp(q0-a/kT) 

[ V M ln{ q j bll }] 

1+exp( — gj-G) 


- {2 e 

b i 


s 






(2 .la) 


where the symbols 0^, X , 0 O , e s » q ’ Y D’ Y n 

and D are described in the list of symbols in Appendix A 
and 6 and p Q represent, respectively the thickness of the 
oxide layer and the fixed charge esntity in the oxide layer 
under thermal equilibrium. 




37 


Let us consider the role of p o in determining the 
harrier height 0^. There are four possibilities that can 
arise; these are positive and negative p Q in the case of 
either JT type or P type semiconductor. Let us consider the 
case of positive oxide charge first. Lor an M-I-H semiconductor 
cell a perusal of eqn. ( 2 . 1 a) shows that with positive p Q , 

0g would decrease. Similar considerations show that with 
positive P Q , 0j for the case of M-I-P semiconductor cell 
increases. Thus, since a large value of 0 -g is desirable 
one would like to obtain if possible, positive oxide charge 
in a p type semiconductor and negative oxide charge in II type 
semiconductor. Technological difficulties in silicon, at 
the present perclude the presence of negative oxide charges, 
thus leaving the only case of interest as positive oxide 
charge in P type semiconductor. An analysis of this has 
been done by Pulfrey [ 9 ] . 

In the case that P Q = 0, exp(q0g/kT) >>1 and 
exp {•(q0 B -E & )/k;T }<< 1 , eqn. (2.1a) reduce to the basic 
quadratic equation in 0 B of Cowley and Sse (see eqn. 21, 
p.374 of [ 17 ] ). As is well known, Cowley and Sze have 
further obtained that in the case that a constant is 
given as 

C 1 = { 2 q e s 11^6 2 / } 

is negligible, the barrier height is given by 
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f>3a " °2 (0 M ~ X) * IU0 2 )( -r - (2 - 1b) 

2 H 
q 6 D -1 

•where Og = (1 + — - — ) 

£ i 

It is interesting to note here that the neglect of the term 
containing C-j is in effect equivalent to ignoring the space 
charge in the semiconductor, which means that a biased case 
such that a flat band occurs is being considered. If this 
consideration is introduced in the beginning itself, i.e. 
in the charge neutrality equation, then eqn.(2.1b) falls 
directly without requiring the solution of a quadratic [18]. 
Once 0g is known, the equilibrium potential drop V ( ) 
across the semiconductor, can be easily calculated by 
knowing V n which depends on the doping and the temperature. 

2.2.2 Illuminated Pells 

1/hen the cell is suitably illuminated from the metal 
side and a load is connected, the band diagram of Pig. 2. 2(b) 
ensues. The following points regarding this diagram should 
be noted. 

(i) Under illumination and ’■>ith a load connected, the 
MIS cell develops a potential difference V such that 

v = V s + V i (2.2) 

where (V^-V B ) and (a - V^) are respectively the potential 
drops across the semiconductor and the oxide layer. Here 
V k is the diffusion potential. Expressions for A and 
are developed in Sec. 2,2.3. Knowing , a value of V g 
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can be obtained from eqn. (2.2) fora given value of V. 

(ii) The barrier height in this case gets modified and 

* 

is denoted by 0g where 

0B “ < VV + V T = <VV + V - V 

= 0 3 + q (2.3) 

(iii) Due to photogeneration of carriers in the semi- 
conductor two quasi-Fermi levels, one for electrons (E^) and 
another for holes (E ) have been shown in Pig. 2. 2(b). It has 

j?p 

been assumed that the majority and minority carrier quasi- 
Fermi levels are constant in the depletion region, a justifi- 
cation of this is given later. Here we note that at a inter- 
face the minority carrier quasi -Fermi level E^, has been shown 
to be separated frem the metal Fermi level E__ by an amount 

of qE . The physical necessity of £ is due to the fact that 
ox ox 

its gradient across the insulator layer provides the driving 
force for the minority carriers to cross this layer. Thus one 
would expect that as the thickness of the oxide layer (5) 
increases, 5 would increase if It is required that 
the same number of photogenerated carrier should cross the 
insulating layer. But note that C OJ£ increases, the 
number of minority carriers at the interface increases, 
which in turn increases the number of minority carriers 
at the semiconductor edge of the depletion region. This 
then modifies the boundary condition for the solution of 
the diffusion equation needed to obtain the profile of 
photogenerated carriers. The net effect of this modification 
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is to reduce the short circuit current Jg^« This happens 
because with constant through the depletion region one 
gets Epp(tj) = C ox j which when used in eqn. ( 1 . 11 b) of 
shows the increases with increasing 

subtracts from J^, Jgg decreases with increasing 5 . 


6 . Since J diff 


(iv) let us now consider the position-dependence of 
Ep n ( x ) and Epp(x) in the depletion region. If the transport 
of the majority carriers is governed by thermionic emission 
then Ep^(x) remains constant in the depletion region [19 1 * 

The nature of Ej,p(x) can be obtained by recalling the 
behaviour of eqn. ( 1 . 9 b), ”hich we rewrite in a modified form 
for N type semiconductor. 


exp [ 


, W b) - 


kT 




1 -1 


= Jp (1/n 2 u p)( ire g /2q 2 kT)^ exp( E pp (o)/kT) 

(2.4) 

Note that (o) = q£ <0. It is apparent from eqn. (2.4) 

■i- jp OX 

that as £ , increases, the difference between the Permi level 

OX 7 

at x = 0 and x = b would decrease. Thus E^ (x) becomes 

jjp 

flatter and flatter in the depletion region as 5 Qx increases. 
In Sec. 2.3.1, we would show that the range of our calculations 
is such that talcing ®jp( x ) constant in the depletion region 
becomes valid. 

(v) The potential difference V at the terminals of the cell 
equal to | Ej, g - Ej^ J /q, where Ep g = E^ in th e bulk of the 
semiconductor. In the conventional situation where £ = 0 

U Jv 
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note that | E p n _E j.p | Al = v * Since the split E-^-Ep is 

controlled by photogeneration of carriers, it is obvious 

that a non-zero value of £ entails a loss in available V. 

ox 

Thus £ is a "necessary" evil for in as much as it 

O X DO 

keeps Jqq constant with increasing <5 but entails a loss in 
the developed voltage available at the external terminals. 

In the case that 6 >6 gg such that Jgg starts falling with 6 
the presence of £ becomes completely disadvantageous. 

OX 

2.2.3 Charge Balance Equation 

* 

The relation 0g = 0g+V^ given by eqn. (2.3) exhibits 

X- 

that to know 0-g one has to know 0^ and An expression 

for 0g is given in eqn. (2.1). Here we setforth to obtain 
an expression for VI. This can be obtained by considering 
the charge balance equation in dark. 


+ %X + Q i + ^SC “ 0 
and under illumination 


(2.5a) 


+ 3 OX + Q i + '^SC " 0 


(2.5b) 


Here Qq X , and Qgg are the charges per unit area in 
the case of thermal equilibrium respectively in the metal, 
the -oxide layer, the surface states and the space charge 
region of the semiconductor. Those v/ith asterisk superscripts 
are the corresponding charges under illumination. In the 
following discussion we would consider = p 5 , Q = p 6 , 

■X* 

where p and p are the charge densities in the oxide in 
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thermal equilibrium and under illumination,, respectively, and 
are assumed to be constants. The functional dependence of 
Qf> and Qgg has been given in Appendix 3. 

Here we proceed to obtain an expression for in 
terms of P, P , Q- , Q. , Q or( and Q.,-, , so that when the 
functional dependences of these are determined, a value of 
Vh may be obtained. The first task is the determination 
6f the potential drop A across the oxide layer in thermal 
equilibrium. 



Solving the Poisson’s equation in the oxide, 
: — § — ; one obtains 


A =¥ (o) -7 (-6 ) 


= [ h 5 p + + Qgg ] (2.6) 

Since A is the difference of potentials at two points, only 
one boundary condition e ^ = (Q^ + Qgg) at x = 0 has been 

utilized in deriving eqn. (2.6). Following a procedure similar 
to the one used for the derivation of eqn. (2.6), under 
illumination one obtains the potential drop across the oxide 
( A - V ± ) as: 

A - Y i = + Q* + Qg 0 I (2.7) 

From equations (2.6) and (2.7) one obtains 
- V i " 


e i [ MQ 0X " Q 0X } + ^i _ + ~ 1 ( 2 * 8 ) 
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Once V. is known V can be obtained from eqn.(2.2) for a 
given value of voltage V developed under illumination for 
a given load. 

So far we have considered under illumination. 

In the ca~e that the diode is not illuminated but is 
forward biased by an applied voltage, let us replace V by 
a corresponding v^; the expression for which can be written 
in analogy -with eqn.(2.S) as 


- -IT [ 4(«ox - Q ox ) + - «i> + «lc - aso 5 1 < 2 - 9 ^ 

f f f 

Here Qq^, Qj_ and Qqq are the respective charge density in 
the oxide, surface states and depletion layer under forward 
bias. 


Thus if the expressions for Q^, Q^, W., Q ort , Q. 


» ''SO ’ ''SO 


and are developed, -one can obtain the values of 7 ^ and 

V^. These expressions have been given in Appendix 3. Note 

that here we have taken a constant charge in the oxide layer. 
In the case that this is not true, one would also require the 

w. jT 

expressions for Qq^» an d ^ox’ surn - iar i ze > it is now 

possible to obtain the distribution of a developed voltage 
V in the oxide region (Vi) and in the semiconductor depletion 
region (V ) . With this distribution known, we can now 
proceed to consider the current voltage relation of the 
Schottky (MIS) solar cell. 3efore concluding let us also 
note that the expressions for charges may themselves depend 
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upon the current voltage characteristics in which case an 
iterative analysis would be required. Two cases where this 
iteration is simpler correspond to the short-circuited and 
open-circuited conditions of the diode. Since the analysis 
in these two cases is simpler and provides physical insight 
in the operation of the device, we have considered than in 
detail in the following sections. 

2. 3 I-V Characteristics 

Our interest is in the calcul ation of the illuminated 
I-V characteristics of an MIS solar cell. Particular attention 
has been paid to two specific cases corresponding to the open 
circuited and short circuited conditions of the cell. But before 
considering these, let us consider the mechanism of current 
transport through an MIS solar cell made on an TT type semi- 
conductor. let J^Cx) and J (x) be the majority and minority 
carrier current densities, respectively , flowing through an 
illuminated cell. The majority carrier current density J n 
is governed by the same equation as in electrically biased 
case under dark and has been discussed by Card and Rhoderick 
Cl 9 3 . They obtain: 

J n = A n T “n ex P(“^ 03 / kr ) I exp(qV/kT) -1] (2.10) 

where A n is the Richardson constant, a q = exp[-X CB (<5) 2 51, 

) is ' tlie barrier height for majority carrier across 
the oxide layer as shown in Pig. 2.2, V is the developed 
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voltage and 0^ is the barrier height under illumination. 
Note that in the case of a dark but forward biased diode, 
is replaced by 0^. 

The minority carrier component is given as 


_ J p - j. L - J diff - J R1 - J R2 


( 2 . 11 ) 


The component is the light generated current density and 
is the same as in the case of S3S0 as given in eqn.(l.13). 

Note that contains the sum of the effect of photo-generated 
carriers in the depletion layer and the quasi -neutral region. 
The component is due to the flow of some of the photo- 

generated minority carriers in a direction opposite to the 
flow of as explained in Sec. 1.2.2. There the expression 
for is given by eqn . (1.11b) which we rewrite for 

convenience: 


J diff " J diff 1 exp * 
where 


kT 


} -1 ] 


( 2 . 12 ) 


J' __ = qD p _/l 
diff H p nO p 


(2,12a) 


Note that we have used the relation Ep^(b) = £ ^ and have also 
introduced the effect of the developed voltage in this 
expression. 

The component J is due to the recombination at the 
interface layer and an expression for this, using Schockley- 
Read-Hall (SRH) theory [20,21], has been developed by Landsberg 
and Klimpke [12]. The steps of its derivation has been outlined 
in Appendix B and lead to the expression 
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Jt>* = J t, [exp{ q( v- 5 0x )/lcT} -1 ] 

.q&L .q0 B - e 


R1 R 
where 

? kT D H(o) { ex P(^)- A K exp( H ^ e - g .)-A > 

In [ 


(2.13) 


J R = n i 


and 


n o (o)(A-B) 


"q$t ' q0| - E p 

{ exp ( -gf) -B > (exp (-V) " S > 


1 


(2.14a) 
2 


A i _ I Qjo) n(o) + H(o ) p(o) ] | |- r G~( o ) n ( o ) + H(o)p(o) H(o)P 0 (o) : 
B ^ on r - — ' t G(o V ' 

(2.14b) 


2G (o)n Q (o) 


2G(o)n 0 (o) 


G(o)n 0 (o)' 


where 

G(o) = T® + T 1 n(o) + T 2 p(o) 
H(o) = Tg + T^ r(o) + T^ p(o) 


The reaction constants T®, T^ , Tg, T^ and T^ are defined 
in Appendix B and will be t sken to be the same for all 
interface states [12], n(o) and p(o) are respecively electron 
and hole concentration at x = 0 and expressions for them are 
given later in eqns. (2.17a) and (2.17b). n Q (o) and p (o) 
are relevant carrier concentrations in thermal equilibrium. 


The component J R 2 is due Bo Bhe recombination of 
the photogenerated carriers in the depletion region and its 
expression is well known I 17 J. It varies as exp [ q(V-£Q^)/2kT ] 
as compared to the exp I q(V- £ QX )/kT] variation of J R1 . Also 
the values of preexponential factors in and are such 

that for the range of our calculations >> J^ 2 * & enCe Bhe 

effect of J ^2 b as Been ignored in the following analysis. 
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Before concluding the considerations of the current 
components, it is important to note that the minority carrier 
density given by eqn. (2.11) must also be the minority 
carrier current density that tunnels through the oxide layer, 
i.e. J can also be written as 

XT 


- J P = V 1 -p 


2 -E+q0 B -q? o x q S 0 X 

a-n exp[ -fi - 2^ ] [1_ exp(- w °-± 


) 1 (2.15) 


“ J tunnel 1 exp(- 4 5 0X /M) -U 
where 


(2.15a) 


a = exp [ -y ( <5 )* 6 } (2.15b) 

P VB 

and 

V 90 B + 4 5 ox "Epp-B, (2.150) 

Here x Tra (5 ) is the barrier height across the oxide for the 

V J3 

minority carriers as shown in Pig. 2. 2. In the following 

sections we use these equations of current transport to obtain 

the short circuit current, the open circuit voltage and then 

the 1-7 characteristics of the MIS cells. The values of parameters 
that have been used are given in Appendix 3. 

2.3.1 Variation of Igg "with <5 : 

The experimental features of Ig^U) an ^IS solar 
cell is that it remains constant upto some 6<dg£ and then 
starts;, falling off [1]. In the following, we have attempted 
to explain this feature. 
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'Ihe short circuit current for an MIS cell on N type 

semiconductor is given as I gc = aJ where a is the area of 

2 

the cell and taken as 1 cm , and is the minority carrier 
current density given by eqn . (2.11), which we rewrite for 
convenience 


• J p ~ J L _J diff" J R1 " J R2 


( 2 . 11 ) 


In the following we have neglected the term due to the 

reasons given earlier. The expression for is the same 

as in the case of Schottky harriers and is given by eqn. (1.1 3). 

2 

We take the value Jt = 25 mA/cm in the present calculations. 
To evaluate and one needs the value of for the 

short circuited case. The general expression for £ can 
be obtained by using eqn. (2.12) and (2.13) in eqn. (2.11) 
and then equating it to eqn. (2.15). The resulting expression 
for ? nT is 


P kT J L +J diff +J R + J tunnel 1 

ox = q < 


(2.16) 


In the case of short circuit, i.e. 7=0, one gets, 


’OX ~ q 


In t 1 +(- 


J diff + J R + J tunnel 


(2.16a) 


Here and are defined by eqns. (2.12a), 

(2.13a) and (2.15a), respectively. 
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We would start the calculations by making two 
simplifications: 

(i) The effect of the component due to recombination 
at the interface is initially neglected. The effect of its 
inclusion is discussed later. 

(ii) The simplified expression for the barrier height 0 B 
in thermal equilibrium, as given by Oowley and Sze [ 16 1, is 
used : 

0 B = O 2 (0 m - x ) + ( 1 -C 2 ) (™ - 0 Q ) (2.1b) 

where 

C 2 = (1 + q 2 6 D s / e i ) -1 

The effect of taking the more general expression for 0 B 
as given by eqn.(2.1a) has also been discussed later. 

We note from eqn.(2.1b) that for the range of 6 

11 2 

considered and for D-jtX'lO per eV per cm one gets 

s 

q0 B - x • have chosen ( 0 -. f ~ x ) = 0.69 eV and 0.5 eV; 

a comraent on this choice is given later. ? ^ is plotted 
from eqn. (2.16a) and is shown in Pig. 2.4(c). In short 

on 

circuited case J (or J on ) represents the total current 

p DU 

SO 

density and has been calculated using the values of Cq-^ 

(i.e. £ nv at short circuit) in eqn.(2. 11). The plot of 

ux 

versus 5 is given in Fig. 2. 3(a) for 0^ - x= 0.5 ©V and 
0.69 eV. Note that starting from 6 = 20 A° a drop in current 
takes place for (0 M ~ x) = 0.5 eV but not for (0 M ~ x )= 0.69 eV. 
The physical reason for this can be seen as follov/s. 



0.50 



A2 69 0 = 



Riots of the open circuit voltage Vo c versus 6. & Plots of the barrier heights, and 
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(i) Let us consider the behaviour of ^ ) 

defined in eqn.(2.12) for V = 0. To start the discussion 

let us calculate the value of ) that is needed io to 

2 

get, say, s ) = 1 mA/cm . For the parameters chosen 

one gets -q £q Y = 0.434 eV. ¥e reiterate that as long as 
-q < 0.434 eV, the value of will not exceed 


diff 


1 mA/cm . Note from Pig. 2. 4(c) that 5^, for (0 M - X )= 0.69eV 


M 


and in the range 0 < 6 < 25 a° , remains below -q C 


OX 


0.434 eV, 


resulting in almost constant value of Jg^Cd). In the case of 

(0j y j- x) = 0.5 eV, -q 5 q-£ reaches a value of 0.434 eV at 

6 = 22 1 ' 5 '/Cgiving J gqff (^sq) = 1 choice of 6g C 

defined above is arbitrary. Here we take it to be that values 

2 

of d where current has failed by 1 mA/cm . Also note that 

2 

in this case at d~ 25A°, = 4 mA/cm , which is comparable 

to the drop experimentally obtained by lillington and Townsend 

[ 11 . 


(ii) The discussion in (i) has exhibited the sensitivity 
of the fall of JqqC 6 ) to the value of 0^, which in the present 
case of low value of D s is given by (0^- x). The value of 
q02 = 0.69 eV corresponds to the workfunction 0^ = 4.7 eV 
for gold vacuum interface and electron affinity x - 4.01 eV 
for silicon. The value ( 0. -- x ) = 0.5 eV has been chosen to 
match the experimental results of barrier height for5>20 A° 
in the case of -&u-Si 02 -nSi diodes [ 143. 
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It is obvious that it is desirable that Jg-< should 

not fall with 6 . The discussion given above points out that 

for this to happen, 0g( d ) should be larger than a certain 

threshold value. This can be obtained by either choosing 

larger value of (0 TT -x ) and/or large value of D ( 5 ) 5 

product. In order to meet the requirement of large barrier 

height at a chosen 6 a large value of 2 0 ( 6 ) is desirable if 

high values of (0,^- x ) can not be obtained. Care, however, 

has to be taken in choosing a large value of 2(6) because 

s 

it also affects the current J-^ due to recombination at 
the interface, which adversely affects the open circuit 
voltage [ 121. Hence, we would return to the choice of 
2 g (6), after discussing the behaviour of the open circuit 
voltage. 


In the beginning we have ignored the effect of 

J Rl ' Let us now include it and see its effect on Jqq( 5 )• 

To include the effect of from eqn.(2.13), the value of 

Ji is needed. Ji is evaluated from eqn. (2.14a). This 
it ri 

needs evaluating n(o), p(o)» n Q (o), p Q (o) ana the values 
of n^, T, G(o) and H(o). The expressions used for n(o) 
and p(o) are 


p(o) = exp [ 


q(0 B - v) 

TT 


n(o) = N 0 exp [ 1 = N n exp [ 1 ( 2 . 17 a) 


‘C 


( 0"R - v ^ 

TT 


kT 
S -E, 


G P S ° X ^ = exp [ 1 (2.17b) 


v 


n 0 (o) and p Q (o) are obtained, respectively, from eqns.(2.17a) 
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and (2.17b) by putting £ cx = 0 and replacing 0 B by 0-g. Note 

that 0g -with illumination but for 7=0 is almost equal to 0g. 

n- is the intrinsic concentration in silicon at I = 300°iC. 

1 

We have assumed G(o) = H(o) following landsberg etal [ 12 ] . 

The expression for G(o) is given in Appendix B. The 


calculation of J R is started by taking = 0. With this 

value of is calculated from eqn. (2.16a). This value 

of *- s now used to calculate J^, and the iteration is 

continued till £ ox converges "within a Unit of error. Using 

the values of 0^ and at 7=0, J^(fi) at short circuit is 

calculated. The maximum value of J^(d ) under short 

—6 2 

circuit is of the order of 10~ anps/cm , hence it is 
justified to ignore its effect on Jgg(d ) since it contains 


2 

Gj/ d ) which has been taken to be 25 mA/cm . Also note 
that when the more general expression for 0 3 given by 
eqn. (2. la) was used, instead of eqn. (2.1h), no significant 
change in Jgg( <5 ) was found. 


While discussing the behaviour of the energy band 
profile at the beginning of this chapter, we had assumed 
that the minority carrier quasi -Perrai level Bp^(x) would 
remain constant throughout the depletion region. We now 
set out to seek the justification for this assumption, 
let us rewrite eqn. (2.4) in a modified form as 


’"®Pp(b ) = “ q ^ QX + kT In t 1+M exp(q £ QX /kT) 1 


Pp 

where M = J N-^ 

hr 


n 


1 

Tv 


)(- 


IT £ 


2q 2 kT 


■)? 


(2.4) 
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We consider the following cases: 

(i) If M exp(q £ QX /kT) <<1, then E pr (b) = £ p , r ; thus 


J Ep'~' “ ^OX’ 


E-o (x) is flat in the depletion region. To get the limit on 

j?p 


1 OX 


for this to he true, we can use the relation 
kT 


^OX = where k <<1. Taking k = 1 and M calculated 

for the parameters mentioned above, we get a typical value 
of q?Q X = -0.036 eV, above which E^ becomes flat. 

(ii) If M exp(q C 0X /kT) » 1, the Ep (b) = -kT ln(M), 
which is the case of S3S0 (where 5 =0). There we have found 
EppCb) = -0.33 el; w.iich is the worst case of the bending 
of Ejp(x). We have already noted an important fact that 
for <5 ) to equal 1 ma/cm , tne amount of Epp(b) needed 

is -0.434eV. This leads us to the conclusion that in this 
case though E^^(x) is not constant, the worst case value 
is such that it does not affect our calculation. 


2.3.2 The Open Circuit Voltage (V^^): 

In this particular case, the sum of the majority and 
minority carrier components of current densities has to be 
zero, i.e. 

J n (x) + J p (x) =0 (2 .18) 


In order to use eqn.(2.18) to obtain I^q one needs to 
get an expression for £q X * This is given by eqn.(2.16) with 
V = V 0 o* Wote that in eqn. (2.16) depends on S QX , hence 
the value of £ QX can only be obtained by iteration. The 
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next task; is to obtain V_.-, . This is done by using eqn. (2. IS) 

(JO 

with J(x) taken from (2.10) and <T(x) taken from (2.15) and 
n p 

K q X taken from (2.16). This results in a quadratic equation 

in exp(qV_. r /kT) given as follows: * 

-E + 2q0 ** 

J L ctp exp( — A 

exp(2qV . /kT) + U-, exp(qV /kT) P 

(2.19) 


= 0 


where 


u i - 


tunnel 

^ J diff +J P 


«„ x&> 4*4 


a h ** 
n n 


( 2 . 20 ) 


Hence one obtains 


1 1, 


In [ - — + ^{ y -j + 


4J t a p expO 


E G +2q0 B ** 


(J diff + J i } 


) K 1/2 

r P > 1 


( 2.21 ) 


Before using eqn . (2.21) for determining Vqq, let us digress 
to simplify this expression and compare the simplified 
expression with the results existing in the published 
literature. 

Let us consider the case where the currents J,. „„ 

diff 

an d ^ are negligible. In this case eqn. (2.19) simplifies 


ezpUVoo) [ J tu nE ,el /M)! 


Ilf* ex P ( 
a n XP ^ kT 


( 2 . 22 ) 
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or 



* 





“L 

7 " 

a 


) 


n 


(2.23) 


Using 0 3 = and the relation n = Vqq/( v oc“^i ^ ’ which 

is the value of the ideality factor for the open circuit 
case; one can rev/rite eqn. (2.23) as 


7 0C - “*b + ^ ln < 


J L • 

A** y a ' 
n n 


(2.24) 


If « n is taken to be unity, eqn. (2.24) is the same equation 
as the one given by Anderson etall 22]which they have 
correlated with several expressions of available in 
literature at that time. Before proceeding further, let us 
consider the nature of variation of 7 00 with 6 as prediceed 
by eqn. (2.23). Here 0 3 and a n are functions of 6 . As 6 
increases a n moaotonically decreases, thereby the second 
term on right side of eqn. ( 2 . 23 ) monotonically increases. 

If one now demands that 7^ fall with 6 , in any range of <5 


it becomes necessary that 0 3 must decrease with Sin that 

range; the- rate of change being of proper amplitude to 

annul the effect of the increase in the second term. We 
* 


know that 0 3 = 0 3 +V j_> and one possibility for decrease of 

* 

0 3 with 6 is that botn 0-g and V decrease with 6 . Perusal 

of eqn. (2.1a) and (2.3) for 0 3 (for (0^-x) < (E & /q-0 o ^and 
(for Q so ?=g^,)--, respectively leads to the following conclusions. 

(A) If (0jyp X ) is constant with 5 then both 0 3 and 7^ 

would decrease if the product I) 6 decreases with increasing <$. 

s 
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This leads to a very important conclusion, that D g with 
constant value, under this model, cannot explain the experi- 
mental feature that decreases with increase of 6 for 
6 > rather it is incumbent that a D g (6 ) function be 

chosen which decreases faster than linearly, for 6 > &q0* 

It is also important to note that to explain the experimen- 
tally observed initial increase of V oc with 6 , the product 
D ( 6)6 must increase in this range. These then are the 

o 

guide lines for choosing a proper variation of D ( 5 ). The 

experimentally observed fact that D g is lower for Si-SiC >2 

system than for bare silicon 1 23 3 is conductive to physical 

existance of such D_ (5 ). The only problem is in choosing the 

proper functional variation, since to the best of our 

knowledge no experimental results in this direction are 

available. To get a physical feel • we have decided to 

2 

study two functional dependences (i) D = a 6 +b and 

s 

(ii) D g = 0 exp(-<$/ g )+d which, as is obvious, have two 

extreme natures of variation of D ( 6 )s . For the sake 

s 

of comparison the results for a constant value of D g 
112 

(= 10 /eV/cm ) have also been calculated and are given in 
Fig. 2. 4(a) . 

(B) So far we have considered (0,,- to be constant with 
6 and required D g to be a function of 6 . The values of 0g 
and V i would also decrease if D g = 0, but (0^- x) tie a function 
of 6. However, this model would require that ( 0^- x) initially 
increases with 6 and th-en starts decreasing, which seems rather 
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contrary to the physical situation. In an actual situation, 
both (0-T- x) a nd D can become monotonically decreasing 
function of 6 , such that the initial increase of 7^ with 6 
is governed by the increase in D (6) 5 and constancy of 
(0jyr~x )» ’•.'hereas tie fall of Vqq for S >6 n( j is due to a 
combined effect of fall of D g ( 6)6 and (0,^- x)* However, for 
the sake of simplicity we have only concentrated on the case 
where (0^- X) is constant and D g (6 ) <5 varies in the fashion 
mentioned earlier. 


To obtain a plot of V ^ versus 6 f using eqn. (2. 21) 

* 

one needs to know J^, J^, J tunne ]_ and 0 . For this an 
iterative process has been used using the following steps: 

(i) and have the same values as in the case of 

the calculation of Ig 0 in Sec. 2.3.1. 

(ii) Taking 0 B and at V = 0, as calculated in 
Sec. 2.3.1, and J^ anne p are calculated. Using these 
values in eqn. (2.21) a value of is obtained. 

(iii) Using this value of for V, 77 is calculated 
from eqn. (2.3) and the relevant expressions of Appendix B. 
Thus 0^ for this value of is obtained. 

(iv) n(o) and p(o) are calcul ted from eqns.(2.17a) and 

* 

(2.17b) for this new 0g. 

(v) 5 at this V = is obtained from eqn. (2.16). 

Then step (iv) is repeated for this new ^ 0 X 
converges within given accuracy. 


until p(o) 
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(vi) and J]; Unae p are calculated for 0g and K q-^ 
obtained from step(iii) and (v) and used in eqn. (2.22) to 
obtain Vqq. Steps (iii) to (v) are repeated until desired 
accuracy is obtained. 

The plots of VqO versus 6 using this proce dure are given 
in Pigs. 2. 3(a) and 2. 4(a). Let us note the following features. 


(i) first of all let us consider the cases for (0-,j~x ) 

= 0.69 eV with D = a <5 ^+b, I) = 0 exp(-<S/ 6 f )+d and 

S o 

11 2 

= 10 /eV/cm which have been shown in Pig. 2.3a(i), 

Pig. 2.4a(i) and (iii), respectively. Note that whereas in 
the first two cases falls off at large values of 5 , in 
the third case there is no such fall off as would be antici- 
pated from the discussion following eqn. (2. 24). Also note 
that these features are supported by the respective plots 
of versus 6 , given in Pig. 2.3(b) and Pig. 2.4(b). 


(ii) The values of Vq C have also been calculated for 
(0pf"X) = 0.5 eV and different distributions of Ds( 5). The 
features of the results are very similar to the previous 
case (e.g. see 2 . 3 ( ii ) ) ; thus from the consideration of 


Yqq alone one may not realize the importance of choosing 
the proper valu e of C0jyr— X ) - However, if the choice of 
parameters is to be such that the experimental nature of 
both Igg( 6 ) and ^qq( 5 ) is simultaneously explained, as 
they must be, there is no choice but to discard the value 
of (0jyj- X ) = 0.69 eV. 



61 


(iii) Regarding the choice between the variation 
2 

D == a 6 + b and D = 0 exp(-6 ,/6*)+ b, one notes that 

from Pig. 2.3a(i) and Pig.2.4a(i), that the variation of 
Y 00 ( 6 ) in the previous case is more atein to the experi- 
mental results [1] than the latter one. This suggests, 
based on the experimental results of Lillington and Townsend 

[ 1 ] that there is a great likelihood that X) ( 6 ) would be 

s 

nearer to a parabolic variation in this case. We have also 
calculated Igg( 5 ) with different functional dependences of 
I>( 6 ) and have found that there i s no perceptible difference 

o 

between the different cases. Hence while discussing )» 

only the case with constant D„ was mentioned. 

s 


2.3.3 Illuminated I-V Characteristics: 

So far we have, concentrated on two particular cases, 
the short circuited ‘.and open circuited conditions, of the 
illuminated I-V characteristics of an MIS cell. But from 
the point of view of efficiency a knowledge of the complete 
I-V characteristics is needed. The basic equations needed 
for this calculation have already been developed in the 
earlier sections; as a matter of fact these have been used 
in the calculation of Iqq and While calculating the 

complete I-V characteristics, the difference from the 
calculation of Ig C and appears in the form of more 
elaborate iterative procedure. This is given below. 


i) Consider a given V at which I is desired. Assume £ 
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* 

ii) Galculate 0^ = 0^+V^ ,by obtaining 0^ ^ rora eqn. (2. la) 

and V. from eqn. (2.3). Since V = V. + V , V is now known. 

X X s s 

iii) Galculate J_ from toe method of Chapter 1 and 
J R1 and J tunnel using eqns. (2.12 ), (2.13) and (2.14). Thus 
the minority carrier component of current is known. 

iv) Calculate from eqn. (2. 16). 

v) Repeat steps (ii) to (iv) till convergence within 
required limit of error is obtained. 

vi) Obtain the minority current from step (iii) and the 
majority carrier one from eqn. (2.10). -Thus the total current 
is known for the given V. 

vii) At each V repeat the steps from (i) to (vi). 


Also, note that if £ ^ is taken to be aero, then, 
obviously no iteration is needed for and this is the 




case treated by Landsberg and Klimpke [12] for computing 
the efficiency in an MIS solar cell. In the case that there 
is interest in the efficiency for the non-zero value of 
£q the iterative procedure given here can be pursued, but 
we have not carried out this calculation. 


2.4 Summary : 

In this chapter an analysis of Schottky (MIS) solar 
cell has been given. She effect of charges due to interface 
states and the effect of the splitting of the minority carrier 
quasi -Fermi level from the metal Fermi level have been *. ■ 
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included in the model. This has been utilized to explain 
the experimentally observed features of an MIS cell (talcing 
Au-SiC>2-nSi cell as example) that the short circuit current 
remains constant upto some 5 and then starts falling, and 
that the open circuit voltage initially increases with 6 , 
reaches a maximum, flattens out and then starts dropping. 

It is worthwhile to compare the present analysis with that 
of Buxo etal t 133. ‘They have found that in the range of 
where Ig^ is falling with 5 , remains constant (see 
Pig. 3 and 4 of [ 131). This result seems to be unphysical 

-Jf 

because this can only happen if 0^ increases with S at 

a required rate; but this is the region where Buxo etal [ 131 

* 

have, in essence, assumed 0^ to be independent of 5 . This 
discrepancy has been removed in the present chapter. In the 
next chapter we undertake to study an alternate method for 
increasing the open-circuit voltage of a Schottky-barrier- 


based solar cell. 



CHAPTER 3 


SGHOTKY BARRIERS WITH OPPOSITELY DOPED INTERFACE LAYER 

A method for increasing the open-circuit voltage of 
a Schottky harrier solar cell (S3SC) is to introduce a thin 
insulating layer between the metal and the semiconductor 
resulting in an MIS solar cell. This has already been 
described in the previous chapter. In this chapter we 
consider another recent technique which increases the open- 
circuit voltage of a S3SC , but as yet has not found wide 
prevalence in the solar cell literature. Due to this reason 
we would elaborate the basic principles of the method in some 
detail. The method consists of oppositely doping a thin 
layer of the semiconductor substrate before depositing the 
metal, thus giving rise to a metal-1 or metal-N-P structure. 
Of course, one can extend the idea to jointly incorporate 
this thin oppositely doped layer along with a thin -insulating 
layer resulting in a metal-(thin) insulator-?-!! or metal(thin) 
insulator-N-rP structure. A perusal of the published literature 
shows that so far there has been only one experimental attempt 
in this direction by Pai et al [ 1 ] , where they have obtained 
the oppositely doped layer by ion implantation and then grown 
thin SiO x before evaporating a thin layer of metal. This 
cell provides some unique advantages which have been described 
later. 
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Before proceeding further, it is important to note 
that the usefulness of the M-I-JT or H-IT-P structures, which 
arises due to their ability to increase the effective barrier 
height, is not limited to solar cells only, as a matter of 
fact, they form a part of the class where an effective 
control, either decrease or increase of the barrier height 
is desired, .fence, in this chapter, we would keep the 
initial discussions at a general level so that these are 
applicable to many devices and later return to solar cells. 

Schottky barriers are employed in a variety of 
semiconductor devices (for example, Schottky-gate field 
effect transistor 121 , Schottky -barri* r solar cells [3,4], 
infection controlled transferred electron devices [41 etc.) 
and in all of these it is desirable to control the height 
of the barrier to optimize the device performance. However, 
it is well known that in commonly employed semiconductors 
like Si, GaAs etc., this control becomes difficult because 
the barrier height cannot be appreciably varied by choosing 
metals with different work functions; rather it is largely 
controlled by the interface surface-state-density which is 
fixed. To obviate this limitation, several workers [ 5-10 ] 
have experimentally as well as theoretically demonstrated 
that the presence of a thin doped semiconducting layer at 
the interface can be used to control the effective barrier 
height. If the surface layer is of the same type as the 
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bulk semiconductor (e.g. metal-N + ( thin )-IT structure) then 
a decrease in effective barrier height is obtained [ 5 ,S J. 

A limiting case of this is obtained when the doping and 
thickness of interfacial layer are large so that ohnic 
contact is formed 15 1. On the other hand, for structures 
like mf tal-( thin)f-N where the interface layer is oppositely 
doped to the bulk, an increase in the barrier height ensues. 
Shannon has obtained the thin interfacial layer, for both 
these cases by ion-implantation [ 6,7} and has compared the 
experimental and theoretical results [ 8]. Recently, 
Van-der-Ziel [ 9 ] has theoretically calculated the barrier 
height for epitaxial metal -I -If structure. The purpose of 
the analysis in the present chapter is to obtain the potential 
profile and thereby the range of geometrical and material 
parameters which would provide the control of the barrier 
height and then apply these to the case of solar cells. The 
calculation includes the effect of mobile carriers which 
allows the determination of some of the paramet c rs which 
would not have been possible otherwise. But before doing 
this, in the next section we give a physical description of 
the role of the interface layer in changing the potential 
profile, briefly review the work of Van-der-Ziel [9};and 
then give the potential profile calculation including the 
effect of mobile carriers. 
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3.1 Pote ntial Profile of Doped Interface-Layer S tr ucture 

let us consider a metal-semiconductor Schottky diode 
which is characterized hy the harrier height 0g o and the type 
and doping of the hulk semiconductor. In the following we 
take the hulk semiconductor to he of N type with a doping 
density The calculations can he easily modified for 

P-type hulk semiconductor. When an interface-layer (P or IT 
type) is metallurgically interposed between the metal and the 
hulk semiconductor, the potential profile ^ (x) in the 
vicinity of the metal contact is modified as shown in Pig. 3.1'. 
If the interfacial layer is of opposite type to the hulk and 
its doping and width are suitably chosen (as discussed in 
Sec. 3.2), a minimum in potential appears at a distance 

x m from the metal (Pig. 3.1h). Thus an effective increase 
in the harrier-height is obtained. If the width of the 
interfacial layer, keeping the doping constant, is further 
increased, there is a flattening of the potential profile at 
the minimum point and the structure behaves as a metal -P- 
semiconductor Schottky diode in series with a P-IT junction 
(Fig. 3.1c). - It is also possible that the width and doping 
of the interfacial layer is such that a minimum in potential 
profile does not appear, i.e. nu harrier height increase 
takes place, as is the case in Fig. 3.1(d). When the inter- 
facial layer is of the same type as the hulk, hut is more 
heavily doped, then depending upon the width and doping of 



r\ 
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this layer one may obtain an ohmic or an injecting contact. 
This distinction depends on two features, the first bring 
the condition whether the space-charge at the high -low 
interface is wiped out or not, and the second being the 
mechanism of current transport, field-emission or thermionic- 
field emission, at the metal contact as discussed in [ 5 3. When 
the interfacial layer is heavily doped and is sufficiently 
thick: such that potential profile of Pig. 3. 1(g) ensues then 
the contact behaves as ohmic. If keeping the doping constant, 
the width of this layer is reduced to an optimum value, then 
the flattening at the maximum of the potential profile 
vanishes as shown in Pig. 3.1(f). This width of the inter- 
facial layer is the minimum width which must be exceeded to 
form an ohmic contact:. In the case that the potential 
profile is monotonic and the doping of interfacial layer is 
such as to lead to thermionic field-emission, an injecting 
contact is obtained. 

We would calculate the general profiles including 
the effect of mobile carriers in Sec. 3.1.2. In the 
following section we consider an M-P-IT structure and 
follow the treatment of '/an-der-Ziel E 9 3 where the effect 
of mobile carriers has been ignored. 

3.1.1 Potential profile of M-P-N diode in depletion 

approximation 

An M-P-N structure is shown in Pig. 3.1(a). Let 
the width of the p layer be W and W be the width over 

Xr ^ 
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which the space charge region extends into the n-region and 
let x = 0 be at the metal contact. Further, assume that 
an abrupt junction has b*e n made. Then in the depletion 
approximation Poisson' s equation are 


d 2 ¥ 

. _ qIT a 



dx 2 

c 

s o 

for 0 < x < W 

P 

(3.1a) 

d 2 Y 

®d 

for W < x < (¥ +¥ ) 

P p n 


d? 

£ £ 

S 0 

(3.1b) 


where q is the charge o± the electron, e g is the dielectric 

constant of the semiconductor and is the permittivity of 
the free space. 


The initial conditions are z 
(a) ¥ (o) = 0,'F ( V w n ) = 0 m1 - 0n . T (3.2) 

where ** m1 and * n denote the vor’i functions of the metal 

electrode and the n type semiconductor and V is the applied 
voltage. 


(b) 

(c) 


¥ (x > and «7dx are continuous at x = V/ 


djfl 

dx 


A , u 

d n 


£ s £ 0 


qN.w 

— 1_E 


P 

(3.3) 


3>F 

dx 


(¥+¥ ) = 0 
P n / 


(3.4) 


The solution for 0< x <_ ¥ is 

=» )W = (ix 2 -x W ) + 2k. v 


e 

s o 


£ s e o n 


whereas for ¥ < x < ¥ +¥ , 
P - p n 


(3.5) 
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qN, 9 q(’^+!'K) ? 

’ (X) - f 2 (X) “ - ^ 1 W P 

(3.6) 

'The width of the space-charge in the N-region is obtained from 
the second boundary condition of eon. (3- 2a), that is 


qlt, P TL+IK 

0 m1"^n- V = * TT“(V W J “ i q(~ -W 


: s e o 


P n- 


(3.7) 


When W n < N a W , a potential minimum exists in front of 
the metal electrode. The position x^ follows from the condition 
dF/dx = 0 at x = x m , which gives 


qN a , qN d 

- a c ( x-W ) + — — W 

e e P e e n 


= 0 


or 


x = 




N, 

W - -2-W 
P N n 

ex 


(3.3) 


whereas its depth A0^(y) is determined by substituting 
x = in eqn . (3.5). 


Using some assumed values of the relevant parameters, 
Van-der-Ziel [g] has computed the amount of increase of 
barrier height which gives physical results. However, we 
have found that though the computations of [91 are valid 
for the range of parameters given in that paper, great 
caution must be exercised when extending this to a different 
combination of the values of parameters. This is mainly 
because of the reason that in 1 9 ^ only a necessary condition 
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for the increase of barrier height has been taken into 

account, whereas another condition needed for sufficiency 

of this has been ignored. This would be discussed in Sec. 3. 2. 

Further, note that fan-der-Ziel has noted that in order to 

get meaningful results for a given II and N , , the width 

st d 

of the 1 -layer should have a lower bound W (min) and 

P 

an upper bound W (max). But his analysis not amenable 

IT 

to direct calculations of W (min) and $ (max), rather a 

ir xr 

trial and error method has to be adopted. A method which 
includes the effect of mobile carriers and allows the 
direct calculation of ¥^(min) and W^(max) has been developed 
in this thesis and is given in Sec. 3.2. But before giving 
this the general calculations including the effect of mobile 
carriers (but in thermal equilibrium) for both and 

M-P-N are given in the following section. 

3.1.2 Potential profiles of H-P-rT and H-II + -N Structures : 

The potential profile 7 (x) can be obtained by 
solving the Poisson's equation: 

^2 = " [ N(x) + P - n 1 (3.9) 

dx os 

where N(x) is net doping density and n and p are, respectively 
the electron and the hole concentrations. The symbols q, £ s , 
e Q have been described in the previous section. The following 
cases where a solution of eqn. (3.9) can be obtained are of 


our interest 



73 


. (a) The effect of mobile carriers is ignored, i.e. 

p=n=0. This is the case discussed in the previous section. 

The calculations done there has inherent limitations as 
■would become clear from the following discussion. As barrier 
height increases, for a P-type interface the valence band 
approaches the Fermi level and the number of minority carrier 
increases. Thus a limit is reached where the calculations 
done under depletion approximation become meaningless. One 
would be interested in knowing this limit below which the 
depletion approximation can be used. It is obvious that 
this limit can only be obtained if minority carriers are 
included in the calculation and the point where departure 
from the results of depletion approximation occurs, is 
noted. This has been our motivation for including the effect 
of mobile carriers in the calculation given below. To the 
best of our knowledge such an approach has not been taken 
earlier. 

(b) The effect of mobile carriers are included but 
only thermal equilibrium condition is considered. In this 
case a closed-form expression for n and p can be obtained 
either by making Ehrenburg' s [11,12] or Boltzman’s approxi- 
mation to the Fermi -Dirac distribution function. In the 
case that one makes Ehrenburg* s approximation, one gets [11,12] 
E -1 

P = 4N y [4exp(-^ + q T)/kT+1 ] (3.9a) 


and 
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n = 4 Nq [4exp(-^ - q y )/kT+1 ] (3.9b) 

where E & is the bandgap, and and are, respectively, 
the effective density of states of the conduction and 
valance band. These expressions for p and n are valid even 
for degenerate semiconductors in case the excursion of the 
Fermi-level in the bands does not exceed 2kT. Hence these 
equations are to be used if any of the semiconductors are 
degenerately doped. Here the reference potential corresponds 
to the Fermi-level in the metal. 

In the case that Boltzmann approximation is used to 
express p and n, one obtains 

p = n.^ expC-q'P/kl) (3.9c) 

n = n^ exp(q'i , /kT) (3.9d) 

where n^ is the intrinsi c carriers concentration. Here 
the reference of 7 is the point where p = n = n-» The 
solutions of eqn. (3.9) with Ehrenburg’ s approximation are 
given in this section. Those with Boltzmann approximation 
can be easily obtained following similar analysis. 

(c) The effect of mobile carriers are included and 
either reverse biased or forward biased devices are consi- 
dered. In these cases, ideally, the current continuity 
equations and the Poisson’s equation should be simultaneously 
solved. However, since the equations are coupled and nonlinear, 
it becomes extremely difficult to get closed form or simple 
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numerical methods of solution, a simplification can be 
brought in by considering the reverse-biased case where 
negligible current flows. Since in this case the potential 
profile can be obtained by solving Poisson’s equation alone. 
Furthermore, one can simplify the situation by -.considering 
either only p or n and not both; for an example in the case 
of M-P-N diode, in the depletion region only p, and not n, 
need be considered. This point would be further elaborated 
at the end of Section 3.2. In the case that a diode is 
forward biased there is no choice but to solve the coupled 
equations mentioned earlier. The only exception is the 
case of an open circuited solar cell, using this structure, 
since in this case again no current flows. 


let us now solve eqn. (3.9) using (3.9a) and (3.9b) 

for the case of thermal equilibrium. We would take the doping 

to be abrupt such that N(x) = for x < x.. and IT(x) = Ng 

for x > x., where x. denotes the metallurgical junction. 

2 3 

Eqn. (3.9) is normalized such that ip= (q? /kT), x = x/L^ 

where 1^ is the Debye length for doping , i.e. 

2 i 

= [kT e e /( q ) ] 3 and the impurity and carrier 

densities are normalized with respect to . w e now restrict 

our attention exclusively to the potential profiles of the 

fnrm of either Pig- 3.1(b) or 3.1(f). In these cases, to 

the right of x.j the electric field is zero for \p= and to 

the left of x • it is zero at \p= # . Using these values to 
0 m 
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eliminate the constant of integration, we integrate the 

normalized form of eqn. (3.9) for x < X ■? a nd x > x . to obtain 

<3 3 



+ 2 2 [ - 


N 2 

^ “ * m) 


F 

4(jp)ln 


1 


P ^ = ^m } 

p(vt=i|i ) 


- 4(5^)ln 


P( y = -4» m ) 
E( v = -4 ) 




(3.10) 


d^ 

dx 


A N _ E( y »4r- ) 

x > x = ± 22 [ -( * R ) - 4(^)ln £- 

X 3 K N 1 F(V = <|> ) 

F, 


where 


E(p =-ip ) , 

- 4(jSi) In ?-] 4 

^1 E( u =_4> ) 

E 

E(y ) = 1 + 0.25 exp(- -V ) 


(3.11) 


( 3 . 12 ) 


and ^ is obtained by solving N^+p-n = 0 in the normalized 
notation and is given as 

% = ln[ (A 2 + (A* + 4A 1 A 3 ) i )/2A 1 ] (3.13) 

where A^= 4exp(+E^/ 2 kT) [1 + ( 4 F v /F^)] 

A 2 = 16 exp(+E / 2kT) + 1 + (16/x^ )(N V -N C ) 

A 3 = 4 exp(+E^/ 2kT)[ UFq/F-j) - 1] 


The signs in eqn. (3.10) are chosen as follows. Eor metal 
P-F diode (Eig. 3.1b) the negative sign is taken for 
X m < X < X-j and positive for 0 < X < X m * The case for 
metal F + ~F diode (Eig. 3. If) is just opposite. In eqn. (3.11), 
the negative sign is taken for metal -P-F diode and positive one 
for metal -F + -F diode. Using the proper signs and noting that the 
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electric field and the potential are continuous at metallur- 
gical junction X = X - > one obtains the potential t|». at x 

J 3 ,3 

by equating eqns. (3.10) and (3.11) as: 


. N ? 

X 3 = ^ _1 ^ [ m “ + 


4 N F(P=.* J 
(^)ln R 


N 


1 


*(u = * *) 


F(v =- O 

+ (4N 0 /N 1 ) In ] 


nr 


(3.14) 


*(**—♦*> 

Ihe use of eqns. (3.9) to (3.14) allows the determination of 
the potential profile in the following manner: 

(i) If the doping N-j of the bulk semiconductor is known 
then is determined from eqn. (3.13). Hote that is a 
positive quantity for donor impurities. 


(ii) A knowledge of Ng and. along with the value of 
if»£, allows the determination of ^ from eqn. (3.14). Here 

1?2 is a positive quantity for donor impurities and a negative 
quantity for acceptor impurities. Methods for choosing the 
value of 4* in particular cases are given in Sections 3.2 
and 3.3. 

(iii) Starting from the value of 1 !'., the next step is to 

cl 

use the laylor-expansion: 




i+1 




(X 


i+1 


V*i 


i+1 




4 


t! 

I 


(3.15) 


where the primes denote differentiation with respect to x • 
It is to be noted that the expansion eqn. (3.15) also 
determines the optimum step-size at each step to obtain a 
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desired amount of accuracy [13]. This is so "because at 
each step ifM is known from eqn. (3.10) or eqn. (3.11) and 
is known from eqn. (3.9). Thus taking the ratio of the 
third term to the second one on the right of eqn. (3.15) and 
equating it to a fixed value, depending on desired accuracy, 
one can obtain the value of ( -X-) at that step. 

(iv) Now the potential profiles in three different regions 
0 < X < X m > X m < X < X^ a nd X > X i are obtained separately. 
Starting from^^, eqn. (3.15) is used to the left till the 
chosen value of is reached. This gives the potential 
profile forX m < X < X y It is important to note that the 
distance ( X. -X m ), which was unknown, has also been deter- 
mined. Repeating the steps to the right of $ ^ till the 
potential reaches ^ R , i.e. the neutral regions, the potential 
profile and the corresponding width of the space charge are 
obtained for x >X^* 

(v) The solution for the region 0< X < X ffi requires the 
value 4 ) = $-r, at x =0. This can be obtained by knowing 

.DO 

the barrier height 0-g Q without the interfacial layer. For 
the reference of potential chosen here, one obtains (see 
inset of Fig. 3.2): 

*50 = ( V 2 * °- 16) 

Starting from ^ :Bo at x = 0, expansion in eqn. (3.15) is used 
to the right till the potential 4, m is reached. Thus the 
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potential profile for the region 0 < x < X ir » along with the 
unknown distance x m » is determined. 

(vi) A knowledge of (x ^ ~X m ) and x m allows the determi- 
nation of xy i^e metallurgical width of interface layer, 
for given values of and \f> ^ so that a desired value 

°f is obtained . It is to be emphasized that the 
determination of this width is a significant advantage of 
the present method of calculation and is of help in designing 
the interfacial layer as is shown in the following sections. 

So far we have considered the potential profiles 
which have an extremum but there is no flattening of the 
profile there. Next, let us consider the profiles given in 
Figs. 3.1(c) and -3.1(g). lo determine the profiles in these 
cases the metallurgical width of the interface layer, along 
with the values of , Ng and^ 3o , should be known, because 
the present method cannot determine the width of the flat 
region near the extremum. Ihe problem becomes basically 
of finding the potential-profile and the width of space- 
charge layer for Schottky barrier between the metal and the 
interface layer, and those of the high-1 ow (or P-N) junction 
between the interface layer and the bulk semiconductor. All 
of these can be easily obtained by the present method. Sub- 
traction of the sum of widths of space charge layers of 
Schottkyxbarrier and high-low (or P-Il) junction from the 
metallurgical width of the interface layer determines the 
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■width of the flat region at the extremum. Thus the potential 
profiles of the type shown in Pigs. 3.1(c) and 3.1(g) can be 
determined. It should be mentioned here that this method, 
in a different context and with Boltzman’s rather than 
Ehrenb rg* s approximation of Fermi-Dirac integral, has been 
used earlier for high-low (or P-H ) junctions with all the 
impurity atoms ionized 1 1 3 3 and with incomplete ionization 
of impurities [14]. However, in each of these cases the 
doping of high side is large enough that Boltzman approxi- 
mation is not valid and the present method should be used 
to get accurate results. 

The remaining profiles where there is no extremum 
but the interface layer is present cannot be obtained by the 
present method. Phis is because the constant of integration 
in obtaining eqn. (3.10) cannot be determined now. 

In the following sections, use is made of the 
calculations given above to determine the parameters of the 
interface layer. 

3. 2 Barrier Height Increase in Metal-P-N Structures : 

It has already been mentioned that an interfacial 
layer which is oppositely doped to the bulk semiconductor 
increases the barrier height, provided the width and doping 
of this layer are properly chosen. This is shown in the inset 
of Fig. 3.2. The considerations that help the determination 
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of these parameters for a metal— P—N structure are discussed 
in this section. Pirst we would consider the case of thermal 
equilibrium and then the cases with applied bias. 

3.2.1 Thermal equilibrium conditions : 

The aim is to determine the parameters such that a 

profile of the type of Pig. 3.1(b), and not of "type of 

Pig. 3.1(c) where a P-F junction is formed or of type 

Pig. 3.1(c) where no barrier increase is available, is 

obtained. This • constrains the value of $ , which was 

m 

taken as a parameter in Sec. 3.1 to lie between <p < ip (min). 

30 m m 

The fact that $ m has to be less than $ for barrier height 

increase is obvious. The value of (min) is shown in 

m 

Pig. 3.1(c) and is determined only by doping F 0 of the 
P-layer at a given temperature. It can be obtained by 
solving (-N a +p-n) = 0 In the normalized notation and is given 
by eqn. (3.13) with F^ replaced by -F . The fact that 
^ ^(min) is determined only by F & at a given temperature 
leads to an important conclusion, i.e. the value of F_ 

cl 

corresponding to m (min) = 4 ' Bq is the minimum value, F a (min), 
which must be exceeded to obtain an increase in barrier height. 

Theoretical calculations withou t recognising the 
existence of F (min) can lead to misleading conclusions as 

cl 

would be evident from the following numerical example. 
Van-der-Ziel [ 9 3 has obtained the potential profile for a 



S3 


me tal-p-n structure without including the effect of mobile 
carriers. He has concluded that for a barrier height increase 
to take place, the inequality: 


EL ¥ > H ¥ 
a p d n 


(3.17) 


must be satisfied, where ¥^ is the v/idth of the space-charge 

layer on H-side of the metallurgical junction. For example, 

1 3 —3 

silicon this is the case for = 10 cm , 0 So = 0.907 eV, 

a sr 'z 

H = 10 cm and W = 1000 £, where W from [ 9) is 

a p n 

8800 1. Por this case, Van-der-Ziel 1 s calculation predicts 


that a barrier height increase of 2 m eV would be obtained 


and the peak in the potential-energy profile of electron 

oo 

would be located 1 20 A in front of the metal. Actually this 
conclusion is misleading because for this case N a (min), using 
eqns. (3.16) and (3.13) with representing -H a (min), is 
1.29x10 cm - which is larger than the value of F & chosen 
here; hence there is no possibility of obtaining a barrier 
height increase. Thus the inequality H a > N a (min) along 
with the inequality eqn. (3.17) forms the necessary and 
sufficient condition for a barrier height in crease to take 
place. Note that the inequality (3. 17)io useful when the 
effect of mobile carriers is neglected, however, the require- 
ment HI, > F (min) is valid in all the cases. 


let us choose a value of H > H (min). The considera- 

3. aL 

tions that help in determination of N & in a particular case 
are given later. Presently our concern is to determine the 
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range W p (min) < W p < W p (max) of the width of p-layer, for 

a given F a such that an increase of harrier height, hut 

without the formation of a P-N junction, is obtained. This 

can he done hy following the method of Sec. 3.1. The value 

of W (min) is obtained when eqn. (3.9) is solved for$ = , f'- 0 . 
v m ho 

Similarly, the value of W (max) is obtained for t (min) 

P mm 

for given value of F . It is expected that the determination 
of W p (min) and ¥ p (max) would he of help in the proper 
fabrication of these devices. 

The general nature of the barrier height increase 
A 0 = 0-g Q > as a function of W p with the doping F & treated 

as a parameter for silicon at 300°K-.. is shown in Fig. 3.2. 
The following features of this figure should be noted. ' 

(a) For increasing value of L, the maximum attainable 
A 0 increases. However, the largest usable value of F may 
be limited by the solubility limit of the dopant, or by the 
fact that as F„ increases the required value of W gets 
smaller making the fabrication difficult, or that (dA0/dW p ) 
becomes sharper so that more precise control of thickness is 
needed. 

(b) For W p > Wp(max), in each case, a constant value 
of A 0 has been shown. This is the region where P-N junction 
forms which is generally undesirable. In conventional calcu- 
lations f 9l, this saturation cannot be obtained; rather a 
monotonic increase in A 0 would continue with increasing W p . 
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This can "be misleading, in as much as instead of improving 
the performance by increasing the barrier height the per- 
formance would deteriorate due to formation of P-N junction. 

Before concluding this section, it should be mentioned 
that the general analysis given in the previous section which 
has been used to analyse M-P-N structures in the present 
section can also be used to analyse M-N + -JF structures. Here 
there are two cases of interest; one of ohmic contact and 
the other of injecting contact. 

A. Ohmic Contacts : In this case the potential profile is 
that of Pig. 3. 1(g) where the space charge at the high-low 
junction is present and provides space-charge-limited current 
to the bulk semiconductor. Since the metal-N + depletion layer 
is separated by a flat potential from the high-low depletion 
region, the presence of the interface layer does not modify 
the barrier for tunneling, hence the value of doping required 
for field-emission would be the same as given in Ref. 15 - 
Regarding the width of the interface layer, it should be 
noted that it must be larger than a critical width, t . , 
to ensure the separation of the depletion layer of metal -ft + 
layer from that of the high-low junction by a flat profile. 
The width t mdn corresponds to the case of Pig.3.1(e), where 
the flat region has just disappeared, i.e. the value of the 
peak in potential is still at ¥ R corresponding to N d2 - It 
is expected that a knowledge of the value of t . would be 
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of interest during the fabrication of Ohmic contacts, since 
it must be exceeded to provide ohmic behaviour. The calcula- 
tion of t min is analogous to that of W (max) in Sec. 3.2, 
hence is not repeated here. 

B* Injecting Contacts : To ensure injecting contacts, the 
doping and the width of the interface layer are chosen so 
that the effective barrier height is decreased by a controlled 
amount. A systematic study of these contacts is of recent 
origin C5,8 land experimental results to achieve injecting 
contacts, particularly in U-InP by epitaxial growth [t6 ] and 
by controlled diffusion of sulfur [17 ] have been recently 
reported. Unfortunately , th e present method of calculation 
of potential-profile fails in this particular case, which 
corresponds to Pig. 3.1(h), as has already been mentioned in 
Sec. 3.1. 

Before concluding this section it must be mentioned 
that the behaviour of the ohmic and injecting contacts alike 
that of the MPN structure, is critically sensitive to the 
applied voltage, hence the present method of analysis which 
is only valid in thermal-equilibrium is of limited interest 
in these cases. In the next section we give an approximate 
method of analysis to study the MPU structure in non-equilibrium 
situation. Similar analysis, with some changes, can also be 
used for the study of ohmic contacts; but for reasons mentioned 
earlier it cannot lie used for injecting contacts. 



3.2.2 Barrier Height increase in H-P-N structure in 
non-thermal equilibrium 

The harrier height increase given in Pig. 3.2 would 
he only valid when no voltage is applied. In the case that 
the device is forward biased the following situation arises. 

The harrier height would increase with increasing 
forward hias. Due to this the number of holes in the P-layer 
would increase and beyond a point one can no longer assume 
that the P-layer is completely depleted; rather a current 
due to minority current injection starts flowing. Thus the 
major advantage of a Schottky barrier diode of being a majority 
carrier device is lost. This limits its main advantage of 
fast switching in switching applications. In solar cell 
applications this means an increased forward current at a 
given voltage which would lead to deterioration in perfor- 
mance. Thus it is clear that the useful barrier height 
increase would be much less than the P-N junction limit given 
in the case of thermal equilibrium. 

The quantitative information about this cannot be 
obtained by solving the Poisson’s equation alone but current 
continuity also has to be solved. This gets to be an involved 
problem due to the coupling of these equations. A decoupling 
can be obtained in the case that the imrefs and 

are assumed to be constant in the space charge region. The 
constancy of 0 n follows from the assumption of majority 
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carrier current transport by thermionic emission and that 
for 0^ follows from the assumption of minority carrier 
current transport by diffusion model [183. It is important 
to note that since the profiles of 0^ and 0^ are already 
assumed, the problem can be solved only by solving the 
Poisson equation and noting that 0^-0^ = 


The solution of Poisson equation itself can be 
simplified by further approximations that p >>n, which 
remains to be valid excepting at a distance of a few Debye 
lengths in the vicinity of the edge of the space charge 
layer. Thus eqn. (3.9), using Boltzmann approximation can 
be written as 


d l _ _ 


dx 


So £ o 


-[ F(x) + n^ exp(-qY/ kT)J 


(3 .18) 


It is important to note the simplicity of eqn. (3.18) which 
has resulted by taking the imref of minority carriers (holes) 
to be constant throughout the space charge region and coin- 
cident with the metal-Permi level even when the bias is 
applied. The majority carrier (electron) imref is displaced 
from the metal Permi level by qV, but since the relation 
p» n has been used, the applied voltage does not appear 
in eqn. (3.13). The effect of applied voltage is taken by 
modifying the boundary condition in the bulk of the 
semiconductor. 
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Let us now solve eqn. (3. IS) for an MPW diode where 
the doping of P and N layers is taken to he abrupt such that 

N(x) = -N» for x < x. 

3 

= Nj) for x > x. 

3 

where x^ is the metallurgical width of the P-layer taking 
the metal semiconductor interface as the origin, let us also 
normalize eqn. (3.18) such that 

ip = , x = x/i^ 

where - , 



and the dopings are normalized with n^, such that 

N 1 = ~ N A / n i for x < X J 

N 2 = Njy/ n i for X > 

¥e assume a potential profile such that a minima ip^ appears. 
This gives a boundary condition d’>/> / dx = 0 where tp = 

The other boundary condition is obtained in the bulk of the 
semiconductor such that d\{>/dx = 0 where ip = where 

is the diffusion potential in the bulk and is given by 
eqn. (3.21). V is the applied voltage which is taken to be 
positive for forward bias. 
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Integrating eqn. (3. IS) once with \p and using the 
above boundary conditions in the proper regions , one obtains: 

for x <X-j ; 

i(d^/ dx ) 2 = t -JJ-jO - ’I' m ) + exp(-<J> )-oxp(-^ m ) 1 (3.19) 

for x> X ^ ; 

2 

Mdi|>/ dx) =[ ^ + 7) + exp (-*|> ) -exp{ -(’I'^-V)}! 

( 3 . 20 ) 

The diffusion potential in the bulk can be determined from 
the charge neutrality condition as 


K = ln [ { 1+ 0 + 

f 2 


(3.21) 


Using the fact that d^/ d.y and $ have to be continuous 
at X = X i ; also using eqns. (3.19) and (3.20) we get 
^i 4> m + N 2 (V R ) + exp -(if H -V) -exp(-4> J 


V 




( 3 . 22 ) 


Now all the relations for obtaining the potential profile have 
been obtained. The calculations can proceed exactly (see 
steps (i) to (v) after eqn. 3.14) as in the case of thermal 
equilibrium given in Sec. 3.1.2. 


It has already been mentioned that the advantage 
of a barrier height increase using a P-layer in M-P-N 
structure gets nullified if minority carrier current gets 
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dominant. An idea about the range where the effect of 

mobile carrier starts dominating can be easily obtained by 

considering eqn. (3.19). Using thi s equation , we have 

2 

plotted (dijj/dx) at x= 0, versus $ in Pig. 3.3. The 

plot, has been given for three values of iL =-U, , which 

along with other parameters of interest is given in Pig. 3. 3. 

The solid curve is for the ease when mobile carriers have 

been ignored, i.e. the exponential terms in eqn. (3. 19) 

have been neglected. The dotted curves are with mobile 

carriers taken into account and the points where they depart 

from the solid curve demark the values of ^ above which 

m 

the neglect of mobile carriers would lead to erroneous 
results. Also note that the points A(i), A(ii) and A(iii) 
on the axis denote the upper limits on ip m for the 
respective dopings due to the formation of P-N junction. 

3.2.3 Open circuit v oltage of solar cells using MPN 
structure : 

Prom the point of view of solar cells the main 
advantage of M-P( thin)->T Schottky barriers over M-IT Schottky 
barriers is due to the increased open circuit voltage 
ensuing from the increase in barrier height. The open 
circuit voltage is obtained by equating the diode forward 
current to the light generated current and solving for the 
resulting voltage. 

The diode current for M-P-U structure can be 


written as: 
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dominant. An idea about the range where the effect of 

mobile carrier starts dominating can be easily obtained by 

considering eqn. (3.19). Using this equation, we have 

plotted (dip/dx) at \ = 0 , versus ^ in Fig. 3.3. The 

plot has been given for three values of Hj which 

along with other parameters of interest is given in Fig. 3. 3. 

The solid curve is for the ease when mobile carriers have 

been ignored, i.e. the exponential terms in eqn. (3. 19) 

have been neglected. The dotted curves are with mobile 

carriers taken into account and the points where they depart 

from the solid curve demark the values of ^ above which 

m 

the neglect of mobile carriers would lead to erroneous 
results. Also note that the points A(i), A(ii) and A(iii) 
on the axis denote the upper limits on i|j for the 
respective dopings due to the formation of P-U junction. 

3.2.3 Open circuit voltage of solar cells using MPU 
structure : 

From the point of view of solar cells the main 
advantage of M-P(thin)-1T Schottky barriers over M-U Schottky 
barriers is due to the increased open circuit voltage 
ensuing from the increase in barrier height. The open 
circuit voltage is obtained by equating the diode forward 
current to the light generated current and solving for the 
resulting voltage. 

The diode current for M-P-fT structure can be 


written as: 




(dy/dX) at 



3 . 3.3 Deviation from depletion approximation (solid line) with mobile 
carriers (dotted line) : Note the scale change in both the cases 
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p 0 -d + A 0-q(V) 

I = aA** I exp [ -q(- — ) ] [ exp(q7/kf)-1 ] (3.23) 

where 0 Bo is the harrier height of M-N structure and A 0 B (V) 
is the barrier height increase at voltage 7 due to the 
presence of P-layer, Other symbols have their usual meaning. 
We assume that 7>> (kl/q) so that unity can be neglected 
compared to exp(q7/kl). Let us laylor expand A0 (V) around 

-D 

A 0 B (o)» where 0 B (o) is ^e barrier height increase in thermal 
equilibrium due to the P-layer. One gets 

dCA 0^(7) ] 

V + .. . (3.24) 

7=0 

Using equation (3.24) in (3.23) and rearranging one gets 


A 0 B (V) = A 0 B (o) + 




dV 


_ . ** m 2 r 

I = aA I exp L 
where 

1 


u( 0 Bo +A 0 b (°) ) 


n = 


1- 


d( A 0 B (7)) 


■] exp(q7/kl) (3.25) 


(3.26) 


~W 


If the value of 


V=0 
d(A 0-r, ( 7 ) ) 


B 

~W 


7=0 


is small compared to unity 


such that n - 1, one can use A 0 B (o) calculated for thermal 
equilibrium to obtain the open circuit voltage. If I B is 
the light generated current in this case one gets: 




A 0 B (O) 



ln(- 


2 

aA I 


(3.27) 


Using the values of A 0 B (o) from Pig. 3. 2, the values of Y^ 

2 ^ 

has been obtained for 0 Bq = 0.503 7, I B = 25 mA/cm ,a = 1 cm , 



94 


A** = 120 Amp/cm 2 /°K 2 and F = 300 °K.. The plot of Y QQ 

versus ¥ > for - 10**^ and 10^ cm”^ and $T^. =,-10^ pad 
1 6 —3 

'110.' eraT djs given in Pig.. 3.4(a)-. 

Recently Oard (18 '3 has done an experimental study 
where he has done a 'solid-state-epitaxial' deposition of 
the thin E-layer on R-type silicon by heat treating eva- 
porated aluminium layers. The experimental results of 
Card [18] are shown in Pig. 3.4(h), The abscissa is the 
temperature of heat treatment., which in effect determines 
and W . Whereas a quantitative comparison of Pigs. 3.4(a) 
and 3.4(b) is not possible due to lack of knowledge of the 
exact experimental parameters, the striking similarity of 
the plots is noteworthy. 

3.3 Metal-(thin)lnsulator-E-R Structures 

A major disadvantage of the tunneling MIS solar 
cell discussed in the previous chapter is that its open 
circuit voltage varies critically with the thickness of the 
insulating layer. Thus the control of the open circuit voltage 
from the point of view of reproducibility from run to run 
needs stringent control on processing. It will be shown in 
this section that if thin oppositely doped layer is grown 
on the semiconductor substrate before growing the oxide, then 
the advantage of the increase of open circuit voltage is 
maintained but its critical dependence on the thickness of 
the insulating layer is obviated. 




Fig. 3.4 Effect of interfacial layer on Voc . 
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In the following we would consider a metal-(thin) 
insula t or -P-2T structure (MIPF); the considerations, for a 
metal-( thin )insulator-N-P structure (MINT) are analogous 
and can he easily followed. A band diagram of MIPW solar 
cell is shown in fig. 3.5(a). 

It is obvious that there are two ways of analysing 
this type of cell. One can either start with MPN structure 
and study the effect of increasing the oxide layer or one can 
start with an MIS cell and study the effect of th e width and 
doping of the P-layer. We have taken the latter approach, 
i.e. we have modified the equations of th e MIS cells given 
in Chapter 2, but used the same analysis of MPF structure 
as given earlier in this chapter. While considering MPH 
structure, for simplicity, we have ignored the mobile 
carriers. 

Let us consider the open circuit voltage of the 

cell shown in Pig. 3.5 (a). This can be obtained by solving 

for V from J +J = 0, where J and J are the majority and 
n p n p 

minority carrier components of current. The modified 
expression for J n is 

* 

q f 0* 4. A 0 \ 

j = A**T 2 a exp( — 2 — ) [ exp(qV/kT)-1 ] (3.28) 

n n 

Note that this equation is the same as in the case, of MIS 
cell excepting the presence of the termA0-g which has been 
added to 0-g. The value of A 0 B for a given doping and width 
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of the P-layer can be obtained from Van der-Ziel’ s analysis 
given in Section 3.1.1. Regarding the minority carrier 
current J , we have assumed it to be the same as in the 

ir 

case of the MIS cell. The expression for the open-circuit 
voltage of the MIPR cell is now obtained following steps which 
are similar to those in the case of MIS cell. The resulting 


expression 


for V 


00 


is 


U + [U 2 + 4J L a P exp f (-V 2 ^ A ^B )/k!r }1 

T ' » ln[ " 3+ 3+ , 

v 00 " q xn 1 9 J 

% 2 (3.29) 

J-exp t(E : + 2qjS* + qA0*)/kI ] + J^ nQel /( (3-30) 


and all other terms excepting 0^ have been defined in 
Chapter 2. Rote that if th e diffusion and recombination 


components are small such that and J R are negligible, 

then eqn. (3.29) reduces to 


In ( J T /A 


(0£ +A0*) 


(3.31) 


Computations have shown that the values of V q from eqns.(3.29) 
and (3.31) are in very good agreement for low values of 6 . 

The open circuit voltage of the MIPR structure as a 
function of the thickness of the insulating layer with width 
of the P-layer as a parameter is shown in Pig. 3. 5 (° and d). 
The calculations have been done using eqn. (3.19) and the steps 
are given below. 
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Let us consider Fig. 3.4(d) first. 

(a) The curve (i) with = 0.0A° is the same as that 
of the MIS cell given in Fig. 2.3a(ii). The material and 
the geometrical parameters that have been used in the 
calculation of Fig. 2.3a(ii) have been also used in the 
following. 

(b) The values of the doping of the P layer and 

cl 

16 3 i c — z 

of the IT layer are taken to be 10 cm"" 3 and 10 5 cm~ 

respectively. The width of the P layer has been taken 
as a parameter and the values are shown in figure. 

(c) at V = 0 is taken from Fig. 2.3. A 0^ is 
calculated using eqn. (3.5), (3.7) and (3.3) for Y = 0. 

(d) Yqq is calculated using eqn. (3.29). 

(e) Then 0g at V=V oC , from eqns. (2.3) and (2.8) and also 
40g (Fqq) from eqns. (3.5), (3.7) and (3.3) are calculated. 

(f ) Steps (d) is repeated to obtain Vqq for the new 
values of 0|j and 4 0*. 

(g) Steps (e) and (f) are repeated until the desired 
accuracy is obtained. 


Using this method the curves of Fig. 3.5(c) and (d) 
have been obtained. In Fig. 3.5(d) the values of W are 

XT 

low so that the limit of barrier height increase due to 

P-IT function formation is not reached. In these cases one 

observes that Y nrl increases with increasing W for all 
Oo p 

values of 6 . But the initial rise of IT is ' still 
sensitive to the value of 6 . Thus the advantages of 
insensitivity of Tqq to ® has not been achieved. 
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The plots (ii) and (iii) of V^, in Pig. 3.5(c) are 
for larger values of W in which case it is apparant that 
the initial rise of has become less sensitive to <$ 
as compared to the MIS cells; thus providing the technological 
advantages mentioned in the beginning of this section. But 
there is one more important point which should be noted here. 


As larger values of are used, the limit of barrier height 
increase due to P-ST junction formation is reached at lower 
values of 6 (see point A on plot (ii) and point B on plot (iii) 
in Pig. 3.5(c)) thereby nulling th e advantages of barrier 
height increase. A calculation of Vqq beyond this point will 
not be meaningful. The analysis given above provides an 
explanation and elaboration of the experimental observation 
of Pai et al [ 1 ] that V qq becomes insensitive to 5 if MINT 
or (MIPN) structures, instead of MIS structure, are used. Let 
us reiterate that this is only true for low values of 6 . 

Pai et al [ 1 ] have obtained the M layer of MINT structure by 
ion-implantating a group V impurity in P type substrate and 
have formed the metallic contact using titanium. 


3.4 Summary 

In this chapter the role of an oppositely doped 
interface layer in increasing the barrier height of a Schottky 
barrier diode has been considered. This property would be of 
use in almost all the devices where Schottky barriers are 
employed e.g. power diodes, microwave diodes, MesPets, solar 
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cells etc. However, the advantage of the in crease of harrier 
height is only available until the minority carrier effects 
start dominating. These limitations have been discussed. 

Prom the point of view of solar cells, Schottky 
barriers have two disadvantages. One is due to the low open 
circuit voltage . This can be obviated by the use of (MIS) 
solar cells or by using structures like MPN or MIPN. These 
have been subject-matters of the previous and the present 
chapters. The other disadvantage of the Schottky barrier 
solar cells, which is also present in the case of MIS, MPN, 
MIPN" and analogous cells, is due to the presence of a thin 
metallic layer through which light may pass. Thus the 
metallic layers have to be chosen to meet the conflicting 
optical, electrical, mechanical and metallurgical requirements. 
A design of the cell, so that some of these stringent require- 
ments are relaxed would be discussed in the next chapter. 



CHAPTER 4 


LATERAL SOLAR CELLS 

In the previous chapters we have studied the 
Schottky harrier, MIS, MPN and MIPF solar cells. The 
difficulties inherent in them due to the presence of a 
thin metallic layer ( 70 A 0 ) on the top has already been 
mentioned earlier. In particular, deterioration of these 
thin metal films due to temperature cycling in actual 
operation may seriously affect the durability of these 
cells, thereby more than offsetting the advantages offered 
by them. Recently, a modified structure for SBSC has been 
proposed by Green [1] in which the requirement of a thin 
metal film has been obviated by using the structure shown 
in Pig. 4.1. This structure is analogous to the grating 
type of photovoltaic cells proposed by Loferski et.al. [2] 
and analysed by other workers [3]. It consists of several 
unit cells, each one of which has an exposed region of 
semiconductor of width 2L and a metal strip of width 
2B = (m-1 )1. Eor brevity, the regions of width 2L and 2B 
are called the lighted and the collecting regions respec- 
tively. Carriers are photogenerated in the lighted region 
and are laterally collected by the Schottky barrier in the 
collecting region. Hence this structure is designated as 
lateral SBSG. The intent of the present chapter is to 
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fig ,4,1 The schematic of lateral SBSC is shown. The width 
of the metal strip is 2B « L(es-I) and the extent of 
lighted region is 2L, The cell consists of sewer a J 
units like this where • the metallic strips are 
connected together. 
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calculate analytically the short circuit current by solving 
a two-dimensional continuity equation and thereby obtain 
the efficiency of this cell. The calculated results .show 
some interesting features. One of them is that the efficiency 
does not increase monotonically with incrase in L/B as 
reported earlier [1 ]. As a matter of fact, it has been 
found that B/B is not a meaningful parameter for this cell, 
rather both 1 and B have to be separately specified for 
uniquely specifying the efficiency of the cell. 

The calculation of efficiency proceeds in the following 
manner. The two dimensional continuity equation for the minority 
carriers in both the lighted and the collecting regions is 
first solved. This is done in the following section. These 
solutions are then used to obtain the spectral response, the 
short circuit current and the conversion efficiency of the 
cell. 

4 . 1 Con tinuity Equation 

In this section the profile of photogenerated 
minority carriers for the lateral S3SC, shown in Big. 4.1 
is determined. To do this let us consider the mechanism of 
collection in a lateral SBSO, bringing out the assumptions 
under which the present calculation has been done. 

(i) Carriers are photogenerated in the lighted region. 

In general, these are transported by drift and diffusion in 
the presence of recombination to the edge of the space charge 
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region (below the metal) which is formed due to the Schottky- 

barrier. Jiowevrr, if it is assumed [4] that the electrostatic 

field is confined entirely to the narrow space-charge region, 

even in the presence of photogenerated carriers, then the 

rate-controlling process in current flow is the diffusion of 

minority carriers. In order to ensure this field-free diffusion, 

the non-equilibrium concentration of minority carriers must 

always be much smaller than the equilibrium concentration of 

majority carriers. This low-level generation needed to neglect 

the drift component of transport, is obtained in the case of 

c-or> 

irradiation of silicon by uncentrated solar light. We confine 

A 


our attention to this ease. 

(ii) The carriers which reach the edge of the space- 
charge layer (below the metal) are transported to the metal 
by drift. Some of the carriers may also recombine in this region, 
We note, however, that [ 4 1 the width (W on ) of this space charge 
region is negligibly small compared with the diffusion length 
(L) of minority carriers and with the width (W) of the cell. 
Typically, in a silicon solar cell Wgg is of the order of a j 

few.microns whereas 1 is of the order of a hundred microns 
and W is of the order of two hundred microns. Hence, on the i 

scale of Tig. 4.1, the edge of the space charge region below 


the metal has been taken at x=0 and any fringing of this 


region in the lighted region has been ignored. In other words, 
it has been assumed that all the photogenerated carriers 
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reaching the edge of the space charge layer are collected by 
the metal. 

Under the assumptions given above the profile of 
photogenerated minority carriers (electrons), in the steady 
state, can be obtained by solving the two-dimensional conti- 
nuity equations 
2 2 

, 3 n. a n. n. 

D ( ryJL + ^2.') _ —JL _ F( Xf x ) (4.1) 

3 * 37 l 

where for the lighted region n^. = and 

F(x, X ) = -Q( X ) a ( X ) exp [ - a ( X )x J (4. 2) 

whereas in the collecting region n. = n and P = 0. Here 

J G 

the reflection losses have been neglected and the quantum 
efficiency has been taken to be unity. The boundary conditions 
(BO* s) in both the lighted and the collecting regions are given 
in Table 4.1. The cell has been considered to be short- 
circuited giving rise to the condition n^ = 0 at x = 0 in 
the collecting region. The determination of the BC's in the 
y-direction poses some problems. This is due to the fact that 
at y = 1, -1, ml and -ml the BC's are not known. To obviate 
this, it is assumed that at these values of y one has n.=f(x), 
and the symmetry of the structure is utilized to put some 
restraining conditions given in Table 4.1. Use of these 
conditions allows one to solve eqn. (4.1) by standard method 
[5 } as shown in Appendix C in terms of the coefficients which 
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Table 4.1 

Table 4.1: Boundary conditions. 


Serial 

Number 

Boundary 

lighted Region 

Collecting Region 

(i) 

x=0 

•n 3n l 

D — == Sn r 

n 3 x 1 

fc* 

o 

li 

o 



(-1 < y < 1) 

-ml < y < -1 

(ii) 

x=W 

= ° 

n 0 = ° 



(-1 < y 1 I) 

(-ml £ y £ -1) 

(iii) 

Boundary 
condition in 

n L = f(x) 

n 0 = f (x) 


y-direction 

at y = +1 

at y=~l and -ml 

(iv) 

Restraining 

3 n j 

L = 0 

3n p 

0 - 0 


condition 

ay 

3y 



at y = 0 

at y = -(m+1 )l/2 

(v) 

Restraining 

3n T 

1 > 0 

3n n 

0 >0 

condition 

9 y 

37 



at y = -1 

at y = -1 



3 n_ 

3 n~. 

(vi) 

Restraining 

1 < 0 

°< 0 


condition 

ay 

sy 



at y = 1 

at y = -ml 
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are to be Iteratively determined. The solution in the 
collecting region is obtained as 


n 0 (x,y) = l C q iin(gyX) cosh(P q (y+(m+1 )l/2)/cosh(E q { (m-1 )l/2] ) 


where fi, are obtained from the solution of 6 W = nir , 
q Q. 

2 2-2 

n being an integer and P q = + L n • To determine the 

coefficients 0 1 s in eqn. (4.3), one would need the solution 

Q, 

of eqn. (4.1) in the lighted region. 


The expression that satisfies eqn. (4.1) and the 
boundary conditions (ii), (iv), (v) and (vi) of Table 4.1 
in the lighted region is: 


ni) (x,y) = g 1 (x)+g 2 (x) 


I Bp sin X p(x-¥) [cosh Qpy/coshQpL ] 
P (4.4) 


g-j (x) = A sinh [ (x-W)/!^] 


Q( X )a ( X) 

gp(x) = — — ? rp — I exp(-aW) exp [ -(x-¥)/L ] 

V a - L n ") 


exp ( - a x) ] 


2 2-2 
Q = X~ + 1 
p P n 

where x_. is the root of the equation: 

Jr 


X p ¥ Cot Xp ¥ + (SW/-D n ) = 0 

The expression for A in (x) is obtained by substituting 
eqn. (4.4) in the boundary condition (i). 

A ={ Q( X )a ( X )/ [ D n ( a 2 - L~ 2 ) ] } (R/T) (4.5) 
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where R =[ ( D n /It n )+S ] exp [Wd" 1 -a ) ] - (S +a D n ) 

I = (D n / L n ) °°sh(W/L n ) + S sinhC V\) 

The coefficients of expansion B , in eqn. (4.4) are determined 
by using the boundary condition (iii) in the following way. 

The coefficients of expansion B^ in eqn. (4.4) and G_ 

P 1 

in eqn. (4.3) are determined as follows. At y = -1, one has 

n O=ni ( : =f(x),i. e . 

CO 00 

I C Sin 3 X = g 1 (x) + g 2 (x) - l B sin x dx-W) (4.6) 
q=1 H H p=1 y y 

Multiplying both sides of eqn. (4.6) by sin jS. n x and integra- 
ting between zero and ¥, we get: 

C n . C 2/w) t-z 1+ z 2 (z 3 -z 4 ) - f (B p 6 n sin x p w )/ ( x p “ 6 n 2 )] 

]p I 

( 4 . 7 ) 

where Z.. = A B [ sinh(¥/l )] /( 1 “ 2 + 8 2 ) 

1 n n ' n n . 

Z 2 = Q( X ) a ( X ) B J I D n (^C 2 - I" 2 ) ] 

Z 3 ={ exp [¥(1^ -a ) ] - (-1) n exp(-aW) > /.( B* + L " 2 > 

Z 4 = t 1-(_1 ) n exp(-a¥) 1 /( a 2 +B 2 ) 

Thus equation (4.7) provides a relation between C n and B^s. 

An expression for B. in terms of 0 ' s is obtained by using 

t q 

the relation 

3 nj / 9 y =3n c /ay at y = -I 
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which gives 

OO 00 

JLl B r ,Q r) sinX T) (x-¥)tanh(Q 1) = I C P sin 3 x tanh(P (m-1 )l/2) 
P P P P P q = 1 Jr " 'x M. 


4 (4.8) 


Multiplying both sides of eqn. (4.8.) hy sinx^.(x-¥) and 
integrating between zero and ¥, we get- 


oo 

B + = Z c l (CP 8 n sinX + W /DC.f - 6?] ) tanh [P (m-1 )l/2 1 


j ^ ^ _ _L _ UJLii A ^ »*/ IA ^ H 

Q. 2 ” *| Xr 'a. M M 


(4.9) 


where Zg = -4 X jJ- { (sin2 x ^ ¥-2 X tanh [ Q^.1 ] } 


Equations (4.7) and (4.9) are iteratively solved to 
obtain the coefficients of expansion B and G inaeqns. (4.4 ) 

P vl 

and (4.3) respectively. Eor the range of parameters for 

which the calculations have been done it has been found 

necessary during the numerical calculation to keep five 

terms in B and fourteen terns in C for reasonable accuracy. 
P y. 


Once these coefficients are known, one can obtain 
the short-circuit current. But before doing this, we would 
digress to consider the case when the surface recombination 
velocity in the lighted region at x = 0 is infinite, i.e. , 
n-^ = 0 at x = 0. Boundary condition (i) for the lighted 
region in fable 4.1, for infinite surface recombination 
velocity, changes to n-^ = 0 at x = 0 for ~L < y <_ L. The 
expression that satisfies eqn. (4.1) and the boundary 
conditions is given by 
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«> COSh Pj 

n L (x,y) = g 5 (x) + g 2 (x) - l B sing x — 9 - 


q=1 


where g^(x) = sinh(^p) 

n 


cosh PI 


and A 


1 


Q( A ) a ( A ) 

D ( a 2 - If 2 ) 
n J 


exp(-a ¥)exp(¥/l )— 1 

[ 2 ] 

sinh (¥/I< n ) 


(4.10) 


For infinite surface recombination velocity (i.e. 
for the case n x = 0 at x = 0) the expansion of (3 can be 

li 1 

obtained in closed form by using the boundary conditions: 


n^ = n^ and 


3 n. 


3 y 


3 n 


3 y 


I . 


at y = -L. 


The expression for 0 thus obtained is 

• uL 

(2/¥) [-Z 1+ Z 2 (Z 3 -Z 4 )1 tanh P q L 
(tanh P I + tanh p ^^" . 1 1 ) 


Thus in this case no iteration is required, this simple 
calculation has been used as a check on the trends that 
have been obtained in the general case of finite surface 
recombination velocity in the lighted region at x = 0. 

Before concluding this section, it should be 
mentioned th<a b recently Hu and ®delburg [3] have also solved 
the continuity equation (4.1). They have used a method of 
Fourier expansion and their analysis has been briefly summari- 
zed in Appendix D. There we have also summarized a difference 
equation technique which has been used by loferski et.all.f2] 
to numerically solve eqn. (4.1). 
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4.2 Short Circuit Current 

The short circuit current Ig^ ( X ) of a unit cell is 
calculated using the expression for from eqn. (4.3) as 
follows : 


w ») - ^ Q 


3 n P 

— -\ ay 

3 x x = 0 


(4.12a) 


00 



WV tanh{ P q (m " l)l/2} 


(4.12h) 


Here unit length in z-direction has been taken. Note that 
the symbol q has been used to denote charge of the electron 
as well as an index of summation. The derivative in the in- 
tegral sign of (4.12a) should actually be evaluated at the 
edge of the depletion region in the collecting region rather 
than at x = 0. However, the evaluation at x = 0 has been 
done as a consequence of assumption (ii) of the previous 
section. 


The spectral response of the unit cell as calculated 
from eqn. (4.2) is shown in Tig. 4.2 for 2B «= 12pm and 
21 = 20pm and 200 P m. The curve for 21 = 20 U m is shown 
for two values ( 1 0 and 10^ cm/sec.) of the surface recombina- 
tion velocity. The values of D , l n and ¥ for which the 
calculations have been done are also shown in -Tig, 4.2. It 
is obvious from Tig. 4.2 that the surface recombination 
velocity at the front surface of the lighted region would 
play an important role in the performance of these cells. 
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To obtain the total short circuit current IgQ» of the 
unit cell, eqn. (4.12) is integrated for AMO sunlight using 
the data of Thekaekara t 6]. A plot of I gg versus the width 2B 
of the collecting region is shown in Big. 4.3. Here 1 has 
been taken as a parameter for different plots. The values 
of D , L and W are the same as given in Big. 4.2 and S=10cm/sec, 

In the earlier literature [1] , it has been implicitly 
assumed that I g0 is only dependent on L and not on B. This 
is at variance with the results shown in Big. 4. 3 where I g( -, 
is also dependent on B. This can be physically explained 
as follows. When B is zero, none of the photogenerated 
carriers- are collected and obviously I gg is zero. When B 
is small, for a given L, some but not all the photo generated 
carriers are collected and I g g increases. When B is 
sufficiently large so that all the photogenerated carriers 
can be collected, then any further increase in B does not 
increase IgQ. 


So far we have considered L as fixed and have varied B. 
It is also of interest to consider the variation of I gc with 
fixed B and varying 1. This is shown as an inset in Pig. 4. 3. 
Initially, as I increases, I g0 increase because photogenerated 
carriers become available. Eventually IgQ saturates when Ii 
becomes of the order of or greater than the diffusion length 
of the minority carriers. These considerations of the results 
in Big. 4. 3 and its inset would be of considerable help in 
explaining the behaviour of the conversion efficiency in the 
next section. 
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Fig. 4.3 Plot of th« •hort-eireroit current vttmi th« width 
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4*3 Conversion Efficiency 

The conversion efficiency of the cell has been 
calculated in the conventional manner [ 7 ] as given in 
Chapter 1. The only care to be taken is that the area, 
for unit length in z-direction, for the calculation of 
the reverse saturation current per unit cell is 2B, whereas 
the total area of the unit cell is 2(L+B).- Thus the 
conversion efficiency ( h ) at a load where maximum power 
is delivered is: 

n = (\ p I mp ) per unit cell/(CP in 2(l+B)) (4.13) 

2 

where for AMO solar radiation P . = 135*3 mw/cm . Y is 

m mp 

obtained from the solution of th e transcendental equation 
qV qV I_ rt 

C + -JSE) exp(-SE) _ (1 + SO) (4.14) 

nkT nKX X S 

where the short circuit current Igg has been calculated in 
Section 4.2, and 

I s = (2B) A**T 2 exp(-q0 3 /kT) (4.15) 

is the saturation current, per unit width, in the collecting 

region. The calculations have been done at room temperature 

2 2 

and for P-type silicon where A** = 32 amp/cm /°K .1 is 
the current of lighted SBSC at the voltage Y . The values 
of the barrier height 0 3 and the diode quality factor n 
which are needed for the calculation cannot be precisely 



determined because they depend on the method of fabrication. 
Hence these have been treated as parameters and several 
calculations have been done. A plot of the efficiency versus 
the width of the metal strip 2B for various values of 1 is 
plotted in Pig. 4.4 for 0-g = 1.1 eV and n = 1.0. 


In order to design a lateral SBSC one would like to 
choose the optimum value of 1 and B. In this regard the 
following reatures of Pig. 4. 4 should be noted: 

1. Por a fixed L and with increasing B, the efficiency 
increases, reaches a maximum and then falls off. The maxima 
is sharp for smaller values of 1 and broadens for higher 
value of 1. It was thought earlier til that increasing 
1/B ratio would monotonically lead to higher effieiency. 
Thus for a given 1 one may think that it would be advant- 
ageous to reduce B as much as technologically feasible. 

Pig. 4.4 shows that this is only true in the region past 
the maxima; in the initial portion of the curves the 
behaviour is just the opposite. Thus, for a given I, there 
is a minimum value of B beyond which it should not be 


reduced even if technologically feasible. We consider that 
this is an important results of the present work. The 
general feature of this aspect of Pig. 4. 4 can be understood 
by considering the variation of Igg, Ig and 2(I+B) with 
increasing value of B. When B is zero the conversion 


efficiency is zero because Igg is zero. When B increases all 
the three quantities Igg» Ig ancL 2(L+B) increase. Eqns.(4.13) 
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width of metal 28 * LCMJU t tot several value# of L. 
The plot of the efficiency ver#u« L for a fixed B i# 
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and (4.14) show that an increase in Igg is desirable whereas 
the increase in Ig and 2(L+B) lowers the efficiency. The 
variation of Ig C has been considered in Big. 4.3 and shows 
a saturation, whereas both Ig and 2(L+B) increase linearly 
with B. These relative changes are such that the conversion 
efficiency, for a fixed L, shows a maxima as seen in Big. 4. 4. 


2. It is also of interest to consider the variation of 
efficiency with I. We confine our attention to that region 
of Big. 4.4 where the efficiency decreases with increasing B. 
Bor a fixed B, the efficiency increases with increase in I, 
reaches a maximum and then falls off as shown in the inset 
of Big. 4. 4. The physical reason for this behaviour can be 
understood by considering the nature of variation of Igg 
with I as given in the inset of Big. 4. 3 and by noting that 
the area 2(I+B) of the unit cell increases linearly with 1. 

It should be noted that again there is a region where the 
efficiency decreases with increasing value of I/B, thus for 
a given B there is an upper limit on I for obtaining maximum 
efficiency. Earlier we have seen that for a given I there 
is a lower limit on B for getting the best efficiency. In 
general, Big. 4. 4, provides a guideline for the choice of 1 
and B while designing these cells and it has to be used 
keeping the technological constraints in mind. 

It is worthwhile to compare the efficiency of lateral 
SBSC with those of conventional Schottky barrier solar cells 
discussed in Chapter 1. Before doing this, it should be 
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noted that in the published literature (for a summary see 
Landsberg and Klimpke 18) ) there is a very wide range of 
effeciencies that have been reported for these conventional 
cells. For example, theoretically calculated values of 
efficiency for conventional silicon SBSC range between 5$ 
to 20$. Such a wide variation is due to the choice of barrier 
height, inclusion or neglect of recombination in the calcula- 
tion, and similar other factors. Hence any comparison between 
the efficiencies of two different types of cells should be 
done, as far as possible under similar conditions. This is 
done below, 

let us consider the lateral SBSO shown in Fig. 4.1. 

If the whole cell is covered by a very thin metallic coating 
instead of the gridded structure, then a conventional SBSG 
ensues. In this case the problem is one-dimensional and the 
solutions, talcing recombination into account, have been given 
in Chapter 1. Considering all the relevant parameters to 
be the same as in Fig. 4.4 and assuming complete transmission 
through the thin metal film, the conversion efficiency of the 
conventional SBSC comes out to be approximately 6.5 percent. 

4- comparison of this value with those given in Fig. 4.4 shows 
that in some ranges the efficiency of lateral SBSC is larger 
than that of the conventional cell, The physical reason for 
this is as follows [1 ]. In the lateral SBSC the light avail- 
able for photogeneration of carriers is less because of the 
presence of opaque metallic layers. Also the total 



121 


recombina tion for carriers is high because of the longer 
collection path. These factors lower the efficiency. However, 
the reverse saturation current of the Schottky barrier is 
also lower because of the reduced metallic area and this 
effect enhances the efficiency. Thus it is possible that in 
some range the lateral SBSC has higher efficiency than the 
conventional SBSC. 

Before concluding this section, let us consider two 
more effects which are of interest in the fabrication of 
these devices, (a) As already discussed in Chapter 2, the 
inclusion of a thin interfacial layer in conventional SBSC 
which in effect forms a MOS solar cell, gives rise to 
increased efficiency. If this is done during fabrication of 
lateral SBSC also, it may be possible to obtain experimental 
values of efficiency which are higher than those of Big. 4. 4. 

In this case, the theoretical calculation would have to be 
modified to account for the presence of the interfacial 
layer as shown in Chapter 2. (b) The efficiency of a lateral 

SBSC is dependent on the value of the surface recombination 
velocity (s) at the front face of the lighted region. Care 
should be taken to keep this value as low as possible. 

4.4 Summary 

The conversion efficiency of lateral SBSC has been 
calculated and guidelines to choose the geometrical parameters 
of this cell have been provided. The deleterious effect of 
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the surface recombination velocity has been pointed out.. 
Perhaps the main advantage of these cells as compared to 
the conventional SBSC would not be in providing better 
efficiency but in their greater durability. In conventional 
cell a very thin homogeneous metallic film is required which 
should have good optical properties, a large barrier height 
and should adhere well to the semiconductor. There are all 
the chances that a thin metallic film which has different 
lattice constant and temperature, coefficient of expansion 
than the substrate would catastrophically deteriorate when 
subject to periodic temperature cycling in actual operation. 
In case of lateral SBSO the metallic film can be made quite 
thick, even an intimate contact using silicides can be 
utilized. These would also improve the series resistance. 

In the light of these observations, it is expected that 
lateral SBSO would be more useful than the conventional 
SBSO in terrestrial applications. 
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conclusion 

The photovoltaic properties of metal semiconductor 
(Schottky barrier) and related devices like: 

(i) Metal- ( thin) insulator-semiconductor tunneling MIS diodes, 

(ii) Metal-( thin)P~N semiconductor and metal- ( thin )H-P 
semiconductor devices, 

(iii) Metal-( thin)insulator-(thin)P-H and metal-( thin )insula tor 
-(thin)N-P semiconductor diodes, and 

(iv) Lateral! Schottky barrier solar diodes 

have been studied in the previous chapters. An analysis of 
the short circuit current and the open circuit voltage for 
these cells has been given. Prom the knowledge of ihe 
complete I-V characteristics in the fourth quadrant (the 
fill factor) the efficiencies of the different cells can be 
calculated as has been demonstrated in some cases. 

It has been pointed out that the temperature variation 
of the parameters of a solar cell are of importance from the 
point of view of actual operation of the cell. Therefore, 
the mechanisms that lead to this temperature variation have 
been studied. Attention has been first given to the study of 
temperature variation of the short-circuit current. It is 
commonly found that the short circuit current increases with 
tanperature. Two mechanisms can be responsible for it. 



(a) The energy hand gap in commonly used semiconductors 
decreases with increasing temperature. Due to this decreased 
band gap some more long wavelength photons are absorbed and 
photogenerate electron-hole pairs thereby increasing the 
short circuit current. Note that this excess of photocarriers 
is due to long wavelength photons which get absorbed deep into 
the base of the semiconductor. Hence, this mechanism would 
only be operative if the diffusion length of the minority 
carriers are exceedingly large. 

(b) The absorption coefficient of the semiconductor 
itself is a function of temperature and it increases with 
increasing temperature at all wavelengths in the case of 
silicon. Thus with increasing temperature, there is extra 
photogeneration of carriers at all the depths in the cell 
and this also causes the increase of short circuit current 
with increasing temperature. Since this mechanism is 
operative at all the sections of the cell, it would be 
important even in the case of materials with small diffusion 
lengths. 

In actual situation, both the mechanisms mentioned 
above would contribute to the increase of short circuit 
current with temperature. However, in the case of a typical 
silicon solar cell, we have found that the contribution of 
the mechanism (a) is almost negligible. Since we have con- 
centrated only on cells made of silicon, the above consi- 
deration is valid for all the cells studied in this thesis, 



125 


hence has not been given for each cell separately; rather 
the numerical results for only the Schottky barrier and 
lateral cells (Appendix D) have been given. Once the 
temperature variation of the short circuit current is 
known, that of the open circuit voltage can be computed 
by knowing the temperature variation of the ideality factor 
and the saturation current which depends on the barrier 
height, A review of the published literature has shown 
that the temperature variations of the ideality factor 
and the barrier height are extremely sensitive to the 
method of fabrication of the devices, hence th e numerical 
computations have not been carried any further in th se 

C&S6S# ■ . * 

A consideration of the Schottky barrier solar cells 
has shown that their biggest disadvantage is the low value 
of the open circuit voltage. This is obviated by using either 
the metal( thin)-P-N and analogous devices or tunneling MIS 
diodes. In the latter case one finds that both the open 
circuit voltage and the short circuit current, hence the 
efficiency are controlled by the thickness of the oxide-layer 
and fall when the thickness of the oxide layer exceeds a 
critical value. 

An understanding of the mechanisms that cause these 
falls would suggest some methods to control the critical 
thickness mentioned above thus help in avoiding the low 
efficiency region. Po-r example, it has been proposed that 
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the fall of the open circuit voltage for thicker diodes is 
primarily caused due to the decrease of the density-of-surface 
states with increasing oxide thickness in the case of MIS cells 
made on N type semiconductor. Thus if this fall is to he 
avoided the rate of decrease of the density-of-surface states 
with increasing oxide- thickness should he kept as small as 
possible. Similarly, the mechanism that controls the cons- 
tancy of the short circuit current with oxide thickness is 
the split of imref of the minority carriers from the metal 
jFermi level, and this should he kept below a critical value 
by obtaining larger barrier heights, to keep the current 
constant. 

The other mechanism to increase the open -circuit 
voltage is to interpose a thin oppositely doped semiconducting 
layer between the metal and the base semiconductor. This 
mechanism- is of very general nature and its uses are not 
limited to photovoltaic devices only. If proper design is 
done, this technique provides a method for increasing the 
effective barrier height which is advantageous in most 
devices where Schottky barriers are employed. Some examples 
of such devices are Schottky gate, field effect devices, 

Schottky clamped transistors, etc. Hence, the methods 
developed for the design of increase of the barrier height 
in solar cells can also be used in the design of several 
other devices. lor this purpose, the graph given in the 
text, showing the increase in the barrier height with the width; 



of the interfacial semiconducting layer and with the dopings 
as parameters can be utilized. Furthermore an analogous 
method can also be used by interposing a thin similarly 
doped semiconductor layer between the metal and the base 
semiconductor to obtain an effective decrease of barrier 
height. This aspect is of interest in designing the injecting 
and ohmic contacts. In the case of solar cells, injecting 
contacts are not^interest, but a consideration for the 
understanding of the behaviour of ohmic contacts based on 
the above principle has been given. Lastly the effects of 
both a thin insulating layer and an oppositely doped semi- 
conducting interfacial layer in increasing the barrier height 
have been amalgamated. The advantageous effect of this 
has been to reduce the sensitivity of the open circuit 
voltage on oxide thickness while keeping the value high. 

In all the cells mentioned so far a very thin 
metal layer is required on the illuminated side so that 
the light can go through. The disadvantages of a thin 
metallic layer from the point of view of long-duration- 
operation of a solar cell under atmospheric conditions are 
obvious. It has been shown that a way of obviating this 
disadvantage is to use lateral cells. The analysis of these 
cells has required a two dimensional calculation. The results 
have shown that the efficiency of these cells are dependent 
on the choice of the areas of both the exposed region of the 
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semiconductor and the metallic collecting region. A design 
graph for this purpose has been plotted in the text. 

Prom the summary of the work presented in the 
present thesis, we come to the following conclusions: 

(a) Schottky harrier solar cells are not promising 
because of their low open circuit voltage. 

(b) The potential of tunneling MIS solar cells has 
not been completely realized so far because these are 
complex devices, in as much as their performance depends 
on cumulative behaviour of several parameters which are 
simultaneously operative. There have been developments in 
this area, but as would appear, the role of all the para- 
meters are not yet well understood. Por example, the 
important role played by the variation of the density of 
surface states with oxide thickness has only been pointed 
out in this thesis. This makes it incumbent that different 
methods of oxide growth which give different variation of 
the _density-of-surface states on oxide thickness be 
experimentally investigated and the suitable one chosen. 

(c)Solar cells based on metal (thin)-insulator-(thin)-P 
-N semiconductor and analogous structures seem to be 
promising, but so far very little attention, both theore- 
tically and experimentally, has been paid to these. The 
most advantageous features of these cells is that they keep 
the open circuit voltage high, but decrease the sensitivity 



of the open circuit voltage on the oxide thickness as 
compared to MIS cells. Thus, these cells would he tolerant 
to less critical manufacturing processes than MIS solar 
cells. 

(d) The manufacture of the lateral cells provide an 
advantage in as much as the choice of the metal for 
collecting junction only depends on the electrical and 
mechanical properties and not on the optical properties 
as in the case for all other cells considered above. Thus 
in the case of lateral cells using MIS contacts,, the harrier 
height and the effective ideality factor would remain to he 
the ones in dark which is advantageous in most of the 
cases. 

To conclude, it must he mentioned that none of the 
cells mentioned above can have an ideal maximum efficiency 
larger than the ideal maximum efficiency available in 
P-N junction solar cells. However, as mentioned in the 
beginning of this thesis the processing of P-F junction 
cells is expensive and there is a search on to locate viable 
alternative structures. We have presented some discussions 
in this direction in this thesis. 
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APPENDIX A 
LIST OP SYMBOLS 

Effective Richardson constant (Amp cm"^ °K”^) 

2 

Area of the cell (cm ), 

Width of the depletion layer (cm). 

Diffusion constants of electron and hole, 

-2 -1 

respectively (cm - /sec - ). 

Density of surface states (states/ cm /eV) 

Band gap ( eV ) . 

Thermal energy (ev). 

Diffusion lengths of electron and holes 
respectively (cm). 

Ideality factor. 

Excess minority carrier (electron) 
concentration (cm ). 

Effective density of states in conduction hand 
(cm ^) . 

Effective density. of states in valance 

-3 

hand (cm 1 . 

Acceptor impurity density (cm ). 

Donor impurity density (cm -3 ) . 

Input radiation power density (W/cm ). 

-2 -1 

Photon flux density (cm /sec - ) per unit 
wavelength ( x). 

Charge of the electron (coulomb). 

Surface recombination velocity (cm/ sec.) 
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I 

V (I ) 
mp v mp' 


qV 


n 




V 


oc 


w 

e_ 


e . 
1 


P 

6 

a ( X) 

hi 

X 

A 


9 < 


Absolute temperature (°K). 

Voltage (current) for maximum power output 
(V(amp) ) . 

Energy difference between conduction band 
and Fermi level (eV). 

Energy difference between Eermi level and 
Valance band (eV) « 

Open circuit voltage (V) 

Potential drop across the semiconductor (V) 
Potential drop across the insulator (V) 
Width of the cell (cm) 

Eree space permitivity (Farad/cm) 

Relative permitivity of the insulator 
Relative permitivity of the semiconductor 
Fixed charge in the oxide (Goulimb/cm ). 
Oxide thickness (cm). 

Absorption coefficient (cm ). 

Workfunction of metal (eV) 

Barrier height (V). 

Electron affinity (eV) 

Potential drops across the oxide in thermal 
equilibrium (V). 

Energy level at surface (V) . 
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APPENDIX B 

CHARGE RELATIONS IN SCHOTIKT (MIS) SOLAR CELL 


In Chapter 2, while discussing the characteristics of 
Sehottky (MIS) solar cell the expressions .for the "barrier 
height and for various charges were required. These are 
given here following an analysis of Landsberg and Klimpke 
(Proc. R. Soc. Lond. , A 354, 101-118, 1977). The symbols used 
h ,€f3 are defined in Appendix A. 


B. 1 Charges Per Unit Area 

Surface state charge : Let us first consider the surface 


state charge, in thermal equilibrium. This is given as 

as, 1 ( B. 1 ) 


0 O 

Vto aJS to) 


where D g is the density of surface states, E is the energy 

of a recombination level at the surface and f . is the 

to 

corresponding ^ermi -Dirac distribution function. If the density 
of surface states £ is constant then Q. becomes 

S -L 

1+exp(q0-g/kT) 


Q. = qD kT tin {- 
i s -j 


+exp (q0 B -E G )/kT 


-> - e! 1 


S' 0 


( B.2) 


Under illumination or when a bias is applied the MIS 
diode would be in non-equilibrium condition. In the case of 
low injection, the departure from equilibrium is not large 
and the system is said to be in quasi -equalibrium. A 
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quasiequlibrium distribution is completely specified by three 
quasi fermi levels associated with electrons, holes and the 
recombination level. Even for moderate departure from 
equilibrium the simplicity of quasi- equalibrium is lost and 
one has to consider recombination kinetics. The analysis of 
one-recombination- center model without trapping and for 
low-injection has been carried out in detail by Schockley 
and Read (Phys.Rev. 87, 835, 1952) and by Hall (Phys.Rev.87, 
387, 1952) and is called the SRH theory. 

They have obtained the expression for the probability 
of occupation of the recombination level as 

G(o)n(o) + H(o) p 1 (o) 

f = ( B.3) 

^ G(o) t n(o)+n^ (o) J+ H(o) f p(o)+p.j ( o ) 3 

Here we have used the notation of Landsberg etal (Op.Cit.) 
n(o) and p(o) are, respectively, the electron concentration 
in the conduction band at the interface x=0 and the hole 
concentration in the valence band at the interface x = 0; 
n^(o) and ( o ) are the electron and hole concentrations 
respectively, when the Eermi level is at the trap energy in 
thermodynamic equilibrium. G(o) and H(o) are given by 

G( o ) = T® + T.nCo) + T 2 p(o) 

(■ B.4) 

H(o) = t| + T 3 n(o) + T^p(o) 

where reaction constants T®, T®, T^ and T^ are 

defined in Pig. B.1 and their values are given in Tabl^- ' 
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Fig.B.1 


Definition of reaction constants occurring in eqn.( B.4). 
The surface state charge under illumination Q* is given hy 

Q£ = lJ g B s f t dE t ( B * 5 


In case, the tunneling from interface states to the metal are 
Ignored, f t in eqn. ( B.5) can be obtained from eqn.( B.3). 


The resulting form of Q* is: 


* 

Q i -T— 


H ( o ) p Q ( o ) exp( - 2q0 B /hT) + exp ( -q0 B /^ ) 

(G(o) n(o) +H(o)p(o) )+ G(o)n 0 (o) 

H('o)p o (o)exp (2(E g -q0 B )AT)+exp(q(E g rfl0 3 AT) 

- ( G(o)n(o)+H(o)p(o>G(o)n 0 (o) 


(G(o)n(o)-H(o)p(o)) s ( exp( q0 B /kT) -A) ( exp( ( 10 B - E g ) % } 

(A-B) G( o)n 0 (o) ’ ln (exp(q0 B AT)-B)(exp((q.0 B -E g )/fcT)-A) 


l D o0c 


J 3 g" ' 
( B.6) 
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With 


(S(o)n(o)+H(o)p(o) ) Gr(o)n(o)+H(o)p(o) , 

+[{ } 

2G(o)n (o) 2G-(o)n (o) 


H(o)p o (o) * 

G(o)n 0 (o) 


( -B.7) 


Although the general expression for calculations is given 
by eqn. ( B.6), the lack of availability of experimental values 
of the reaction parameters and their variation with oxide 
thickness in' both the cases of an ’illuminated* and ’biased’ 

MIS diode, makes it difficult to obtain me' ningful quanti- 
tative values. Here for the special case$ of G-(o) = H(o) 
expression for Qt has been developed. The advantage of the 
approximation G-(o) = H(o) is that the unknown recombination 

parameters cancel out. If we further make an approximation 
2 

[ ] » [ p Q (o)/n o (o) ], eqn. (A. B.6) reduces to: 


* 

Q i = 


qB kT 


tin t 


exp(q0 B /kT) + ( p(o)+n(o) )/a Q (o) 
exp((q0 B -E g )/kT) + (p(o)+n(o) )/ n Q (o) 


E 

3- dh 


[ 1 + 


q(o)-p(q) 

n(oT+pTo)' 


3 } + qB 0 


( . B.S) 


It is interesting to note that if p(o) »n(o), then eqn.( B.S) 
simplifies to: 

< - 0 


Let us now consider the interface charge in the 
case that the diode is forward biased, but there is no 
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f 

illumination. The expression for ^ is the same as that 

•Wr 

of in eqn. ( 3*6) excepting that 0^ should be replaced 

f 

"by 0B an< i n(o), p(o) should he accordingly calculated. Also 
note that the reaction constants would be different in the 
two cases. 

Space charge in depletion region : The space charge is the 
depletion layer Q* G is 

<4 = iv - 12 e s <i vw 4 ( b - io > 

1 

by the depletion approximation with w - [ (2 / 2N^) (V -V g )]^. 

In equilibrium V g = 0 and Q SC is given by 

«so = t2 l ur B v b ji ( B.11) 

. B . 2 : Potential d r ops and Barrier Height 

The potential drop V^-V s across the depletion layer 
is given under illumination by 

o p<5 2 <$ 

A- V i = J E i (x) dx = —g-- -- + (Qi+Qgo) — ( B.12) 

Thus , 

(e./6 )( A- V i ) = ^ 5p + Q L + Q s0 ( 3.13) 

becomes in equilibrium, • 


(e i / 6) A = i 6p o + + Q sgq 


( B.14) 
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Subtracting, 


-( ® ± /*)\ = 4 < p - "„) ♦ CQi-Qio) * .COso-Soq) { B ' 15) 

Now£ from 0 B = 0 M - X- A , ( B.11), ( B*2) and 
(B.T3')» the equilibrium barrier height is given by the 
solution of 


0 B - 0M- x ‘ 



2 e. 


l 


J s <1 

e . 
x 


^ e G*1cT In { 


1+exp(q0-r,/kT) 


1+exp(- 


w 


*5 

kT 


} J 


) 


o n 2 

ii- 1 0 

6^ S^O 



e . 


l 


2 e Q q 


wv 


( B.12) 


The values of the parameter used in eqn. ( B..12) are given 

in Table •. B.l. Also note the values of tunnel exponent 
used are given in Table B.2. 
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Table B.1: Data Adopted for MIS 3olar Cell*. 


Parameter 

Numerical 

Dimensions 


value 


1.12 

eV 


4.01 

eV 

T 

300 

K 


3x10 l7 q 

■O' cm - ''* 

t®=t| 

1 . 1 2x1 0~® 

3 -1 
cm s 

Ti=T 2 =T 3 =T 4 

3.66x10' 25 

6 -1 
cm s 

A**= Ap* = A** 

120 

A c :ji~ Z K~ 2 

e 

l 

11.8 e 

0 

P/cm 

8 

S 

3 - 9e o 

Ff cm 


12.5 

-2 -1 

cm sec 


100 


% 

io 15 

cm~^ 

n i 

1 .6x1 0 10 

cm' 3 

* Taken 

from Landsberg and Klimpke (up.CiU. ). 

Table 

B.2: Tunnel Exponent X 2 

6 

«<X) 

5 10 15 

20 25 


0 


2 


4 


8 


14 
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APPENDIX 0 


SOLUTION OE TWO DIMENSIONAL CONTINUITY EQUATION 


0.1 Solution in Dark 

Let the carrier concentration in the collecting region 
he n n (x,y) where 


n n (x,y) = XY 


( 0 . 1 ) 


Let the separation constant he X= + 8 . Then one 


obtains: 


1 d 2 X o 2 

v — 2 = ± ® 

X dx 


( 0 . 2 ) 


1 d Y _ J 2 

Y T~2 2 - ± e 


(0.3) 


It can he easily checked that for X = 8 there is no 

2 

solution for eqn. (0.2). For X= - 8 , we have 


X = t cos 8 x + r 1 sin Bx 


(0.4) 


Y = G exp ( /P x } + D exp { -/P x } 


where 


P = + 6 2 

L h 


(0.5) 


( 0 . 6 ) 


Use of B.C.’s X = 0, at x = 0 and W shows that the solution 
of eqn. (C.4) is 


X = I t’Sin B n x 

q=1 q q 


(0.7) 
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Use of the B.C.’s dn c /dy = 0 at y = - |(m+1 ) shows 
that 

coshj /P{ y + ^(m+l ) } j 
cosh[/P |(m - 1) ] 

where K = D exp(/P §(m+1 ) )cosh{/P |(m-1 ) } 

It has been checked that the complete solution in the 
collecting region, i.e. n^C^y) is obtained as the product 
of (0.7) and (C.8) and is given as eqn. (4.3) of the main 
text. 


0.2 Solution in the Lighted Region 


Let the carrier concentration in the lighted region 
be n L (x,y) and expressed as 


n L (x,y) = X^x) Y-j (y ) (0.9) 

2 

In this case let the separation constant be X = + x .It 

2 

can be shown that for X = x there is no solution. In 

2 

the case X = - x one obtains 


X^(x) = A cos x x + B sin xx (0.10) 

(y ) = E exp (vQ y ) + P exp(- / Q y ) (0.11) 

Use of B.O. X^ = 0 at x = W, gives from eqn. (C.10) 

A = -B tan x^ (0.12) 

dX. 

Use of H, ^ SX 1 at x = 0 gives for B ^ 0 
X cot X¥ = -S/P^ 


(0.13) 
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dY. 

The use of B.C.; = 0 at y = 0 gives 

E = -P (C.14) 

m 

Note that eqn. (G.13) has multiple solutions. Use of this 
fact and eqn. (C.12) and (0.14) in (0.10) and (0.11) gives 
the third term of the solution of nj j (x,y) as given in eqn. (4.4) 
of the main text. 

The terms g^ (x) and g 2 (x) of eqn. (4.4) are obtained 

2 

as follows. When A = 1/Ir tt > one has the solution g^(x) of 
eqn. (4.4). The term g 2 (x) is the particular integral due 
to inhomogenuity of the continuity equation. We have guessed 
the particular integral to be of the form: 

Q( A ) a ( X ) 

g 2 (x) = -~2 r2 { A exp [ (-x+b)/ 1^1 - exp(-<*( x )x) * 

Dh( a _ ) 

( 0 , 15 ) 

Using the boundary conditions, g 2 (x) = 0 at x = W, in 
eqn. (0.15) we obtain the g 2 (x) of eqn. (4.4) of the main 
text. 
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ALTERNATIVE ANALYSES OE IATEP. CELLS 
D. 1 Analysis of Lateral Cell by Fourier Series Expansion : 


In a recent paper [Solid St. Electronics, 20,119,1977), 
Hu and Edelburg have analysed the lateral cell structure of 
using the technique of Fourier series expansion instead of 
"boundary matching" as has been done by us in the previous 
sections. We summarize, in the following, the relevant 
portion of the analysis of the Hu et al so that a comparison 
with the "boundary matching"can be done. 

The differential equation in the two cases is the same, 
as given by eqn. (3.1). However, the boundary conditions that 
have been utilized are different, in as much as, in Hu's 
analysis , the B.C. number (iii) in 'lable 3.1 is not needed. 
Also note that we have analysed the structure based on 
a P-type semiconductor, whereas Hu etal have done it for N-type 
semiconductor. Their solution is written for holes WhiLch are 
the minority carriers. 


The solution is expressed as sum of particular 
solution Pp and the homogeneous solution p^. The solution 
Pp is given by 


p„ =xaQ . „ £_ 
r p m=o 


' a d V G m^“ e " a d ' lcos X 

"2T2 72“ 


m 


1- L 


a 


+ L 


V 


m 


(D.1) 


where L is the diffusion length in centimeters, = ihtt/Cm) 
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U Q = 1 and u m = 0 for m ^ 0, 

v 9 sin (3 X ) 

0 “ b ’ G m ^ L+B ' ~TZ 
m^o m 

oo sh 6 ¥ - 

¥ m = e" ax - cosh 0 x + ( — r )sinh B x 

m . m sin 3 m W m 

f 1 2 v-4- 

■where 8 m = + \ m J 

1 

and- d is the thickness of metal stripes. 

The related homogeneous solution is obtained by 
the technique of eigenfunction expansion. 


P h = I A m sinhC $ m (W-x ) > cos^y /sin 8 m ^ 


(B.2) 


where the A^ 1 s are the Fourier series coefficients of the 
excess hole concentration along the boundary x = 0 and A m l s 
are given by the following expression. 


\ = - " d )!J 


+ (l_u ) ( 0 4 G ) t A R + ( 1 -e~ u ) 
v nr v o 2 m' mm ' 


Q 0 

- °t d\ m m ] 


fi _ i u ) I (0 . + C .) f A.R. + ( 1 -e ) 4 — A 3 

^ 4 u m ' ^ v m+3 m-3' 3 3 b 

dM 

(D. 4 ) 


where R = -(D/S) & coth & ¥ 


1-L‘ a + 1 X' 


cosh B ¥-e 
^m^ sinh B ¥ 


-a ¥ 


■) - a ] (£*4) 
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Note that eqns.(A.5) and (A. 7) in Hu etal’s paper have 
algebraic errors which have been pointed out by us 
and the corrected form is given in eqn. (D.3) and eqn.(D.4). 
Ihe infinite set of equations (H. 3) is modelled as a system 
of m simultaneous equations to be solved for the m lowest 
order Fourier series coefficients A^. At least m = 70 
coefficients were needed by Hu etal for reasonable accuracy. 
Once the solution for the minority carriers are obtained, 
the calculation of efficiency etc. in the two methods is 
similar hence would not be discussed again. Thus so far, 
the comparison of the two methods is concerned it is to be 
only done from the point of view of ease of numerical 
computation in obtaining th e minority carrier profile. Our 
experience shows that solving for several simultaneous 
equations would definitely be more horrendous than the 
computations needed in the "boundary matching" methods. 

Before concluding this section we would also like to 
point out that though the surface recombination velocity 
at the exposed part of the semiconductor has been incorporated 
in the analysis of Hu etal, its physical significance has 
not been explored at all nor have been the, design aspects 
emphasized, keeping an appreciation of the operative physical 


mechanisms in view. 
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D.2 Numerical Analysis of Lateral Cell 

The disadvantages of a dead layer on the top surface 
of a P-N junction solar cell were noticed by many workers 
but perhaps the first method of obviate this, was proposed 
by loferski etal [ Technical Report NGR-40-002-093/2 to the 
National Aeronautics and Space Administration, June 1975^* 
His suggestion was to use a grating structure where "unlike 
a standard diffused cell which has its front surface covered 
completely by a thin diffused skin region, a grating type 
cell would have part of its front surface covered by a 
grating of finely spaced fine grid lines which form a 
potential barrier with the base substrate. The grating 
structure junction could be a diffused P-N junction, an 
alloyed P-N junction or a Schottky barrier junction". Thus 
we see that Loferski’ s grating structure is the one which 
we have called the lateral cell. He has set up the two 
dimensional continuity equation which is the same as our 
eqn.(3.1) and has solved it numerically by using finite 
difference technique. For the solution of the simultaneous 
equations he has used the successive overrelaxation method. 
As is well known, the disadvantage of a numerical method as 
compared to an analytical solution is that each problem has 
to be solved on its own and general design guidelines cannot 
be obtained. Hence we feel that the analytical method given 
in Section 3.1 would be more convenient than the numerical 
method given by loferski. 
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D.3 Lateral Schottky barrier Solar Cell - Simplified Model . 

The analysis of lateral SBSC given in the previous 
section requires computation of the matching coefficients 
Bp* s and C ’ s which, though rigorous, is time consuming. 

In many cases, it would be desirable to develop an approxi- 
mate analysis which is computationally simple. Such an 
analysis of the lateral SBSG is provided in this section. 

The simplification is brought out by the fact that in this 
case no iteration is required while solving for B ' s and 

C «s. 

q 

The method consists of solving the continuity equation 
in the collecting region as done in the previous section and 
of modifying the continuity equation in the lighted region 
as follows: 

,2 

d n T n T „ _ 

3) —T-- — = -Q a e” X (D.5) 

dy" x n 

Here the diffusion component in the x-direction has been 
neglected. Hence the boundary conditions (i) and (ii) in 
the lighted region in Table 3.1 are not required. The solution 
of eqn. (3.16) is obtained as 

n T = 2Au cosh-^— + t Q a e ax (L.6) 

II I n a 

where A.j is the coefficient to be determined later. 

In the collecting region the solution of the two 
dimensional continuity equation is given by eqn. (4.3) as in 


Section 4.1. 
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Using the condition that y = -L, n^ = n^ = f(x) 
and 3^-,/ 3y =3 n^/3 y, the coefficient CMs, and A^ can 
be expressed as 


2t Qa(3 (l-(-1) n e " aW ) > 

^ _ n n 

n p p P (in— 1 )L j 

¥( a +3 ){ 1+PLtanh ( )coth~ } 

n' n n c l 

n 

and 


(D.7) 


A 


1 


2sinh(L/L n ) 


I 0 sin(£ x) P tanhfj 2=1 j,) 
q_1 y h a y ^ 

(D.8) 


Onc-e the coefficients ^ f s are known the expression of short 
circuit current is given by 


> 

(D.9) 

Note that the symbol q has been used to denote charge of 
the electron as well as an index in the summation. The 
particular use is clear from the context. 


I s0 ( X) 


4q Q« L' 


TT 


n 


oo 3 *■{ 1-(-1) n e~ W }tanh{ } 

E -. hi , „ u ' — ' ■■ - jy- , — - 

2 P S 

q_1 P„(« + 3 2 ) (1+1 P _ t anh ( — S- ) tanh y 
q q n q 2 h 


The plot of Ig 0 for AMO sunlight versus the width 
of the metal stripe for several values of 1 is shown in 
Pig. P.1. In the calculation of eqn.(D.9), fourteen terms 
in the summation have been taken. The following points regarding 
the plot of Pig. D.1 should be noted. 
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(1) It is worthwhileto compare the results of the approxi- 
mcite model as given in Pig. D.1-, with those of two dimensional 
calculation given in Pig. 4. 5. One notes that the basic trend 
of the results in the two cases is the same. However, the 
numerical values are different and the amount of variation is 
given later. Here we only want to point to a general trend 
that is always present. This is that the value of Ig 0 found 
from the approximate calculation is always lees than that 
from the two dimensional calculation. The reason is the use 

of one dimensional continuity equation in the lighted region. 

In this case the effective length that the carriers have to 
traverse before collection is increased, thereby increasing 
the chances of recombination thus lowering the short circuit 
current. 

(2) The amount of difference between the graph of Pig. 4. 5 
and Pig. D.1 depends upon the values of 1 and B being consi- 
dered. One finds that for 1 * 50 ym for the whole range of B, 
the difference between the two calculations is less than 20$. 

In case one is only interested in getting an approximate idea 
but with lots less numerical effort it is obvious that the 
approximate method leading to Pig. 4. 3 would be more convenient 
than the two-dimensional analysis. Por exaaple, one such use 
can be envisaged when one is interested in studying the 
temperature variation of the short circuit current as mentioned 
in Chapter 1 . 



155 


Such calculation of variation of with temperature 
using the expression for the absorption coefficient given in 
eqn.(1.2Jl) for silicon, is given in Fig. D.2. 



0.20 



= 32 pm 


o 



Fig= D.2 Temperature variation of short circuit current with 1=300 /am and B as 
parameter. 
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